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On Steady-State Solutions 
of the Navier-Stokes Partial Differential Equations 


ROBERT FINN 


1. Definitions and Summary 
A three dimensional vector field u(#) =u(x,, x), ¥3) defined in a region & of 
Euclidean space will be said to represent in & the velocity field of a time independent 
flow of a viscous incompressible fluid, provided, a) all derivatives of u(a) up to 
the second order exist and are continuous in &, and b) there exists a scalar (a) 
with continuous first partial derivatives, such that the Navier-Stokes partial 
differential equations, 


(1) Au—Vp=(u-V)u, 


Givity— Or 


are satisfied by the functions u(x), #(x) in &. The notation conforms with the 
classical usage of vector analysis. Certain constants which appear in the usual 
formulation of the equations (1) have been arbitrarily set equal to unity, because 
they play no role in the context of this paper. 

Such a vector field u (a) will be said to define a solution of an exterior boundary 
value problem provided a) and b) are satisfied, and provided c) & is the exterior 
of a bounded region G whose boundary surface 2 consists of a finite number 
of components 2}, ..., 4',,, d) on each component 2;, u(x) assumes prescribed 
values u**, and e) u(a) tends to a prescribed constant vector uy as > oo. It 
will be assumed throughout that the surfaces 2; and the prescribed data u* 
are smooth. For precise assumptions, sufficient in the present context, see 
Opevist [/]. We shall also make the (physically natural) assumption on the 
boundary data w;*, that the surface integral 


(2) Jux-dS =0 


for each component 2;. We note in particular that (2) is satisfied when uy =0, 
7=1,...,m. A solution u(a) corresponding to such data admits a physical 
interpretation as a flow of a viscous fluid past a rigid obstacle. 

In a remarkable and penetrating study of the Navier-Stokes equations, 
Lreray [2] has proved that for any region & of the type considered, there exists 
a solution of the exterior boundary value problem in a generalized sense. Precisely, 


This work was prepared while the author was a Guggenheim Fellow, 1958—1959. 
The author is indebted to the Guggenheim Memorial Foundation, and also to the 
California Institute of Technology, for support and encouragement offered in many 
ways. 
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LerRAY has demonstrated the existence in & of a vector field w(a) satisfying 
a), b), c), d), and such that w(a) is, in every compact subregion of &, the uniform 
limit, as No along a suitable sequence of values, of solutions wy (a) of the 
following interior boundary value problem: ty (x) =uj (x) on 2;, Uy (x) =U On 
the boundary Sy of a sphere Sy of radius N and center at the origin of coordinates. 
Further, the Dirichlet integrals of uy (a) over the regions éy=SyM€& are uni- 
formly bounded, 


(3) D[tty] = f |Vuy|24V <K 
én 


where K is independent of N. 


LeraAy did not prove that his solution u(a) tends to Uy as a limit except in 
the special case Uy=0, a case which seems to be without physical interest. The 
chief purpose of the present paper is to fill this gap. We shall prove the existence 
of a function R(¢) such that for any ¢>0, the limit function wu (a) has the property 
| (a) —uy|< whenever |x| > R(e). We shall also obtain an integral representa- 
tion which is valid for any solution of an exterior boundary value problem 
obtained by the method of LeRAyt. Our procedure exploits a reformulation of 
the equations (1) suggested by a linearization due to OSEEN [3], and a form of 
Green’s identity used by Opgvist [J] in his study of the linearization of STOKEs. 
We are led to an integral representation for the solution #, which differs slightly 
from the one employed by LERAy and which proves to be advantageous for 
obtaining the necessary a priovi estimates on the approximating solutions. We 
shall prove, then, the following result. 


Theorem 1. Corresponding to any region & of the type considered, and to any 
prescribed data Uy and ux satisfying (2), there exists at least one solution of the 
exterior boundary value problem a), b), c), d), e). The pressure p (a) corresponding 
to this solution tends to a finite limit py as e-> co. The solution and pressure admit 
the representations 


ct (@) = Uy — f {X(w, y) - [ue(y) — uy] - Vuely)} ay + 
é 
(4) +f {Xw,y)- Tuly) — (ely) — uy] - TX (a, y)} dS, — 


— Jf [w(y) — uy] -X a, y) (9° n)dSy, 


p(x) = po .. {p (x, y) - [w(y) — ug] -Vuly)} avy + 
(4’) as a { (x, y) - Tuy) — [u(y) — uo] - Tp (a, y)} dS, — 
— J [e(y) = uo] (ae, y) (m9) d Sy. 


Here ny denotes a unit vector in the direction of u,. The tensor X and vector p 
are explicitly known and will be defined in the following sections, as will the 
operator T. 
| As a corollary of our method, we are able to obtain a characterization at 
infinity of a class of solutions of (1) which includes the solutions of LERAY. 


1 The representation is v 
cf. Theorem 2. 


alid for a class of solutions which may be more general, 
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Theorem 2, Let u(x) denote a solution of (1) in a neighborhood & of infinity, 
and suppose that the Dirichlet integral of u(x) in € is finite, [|\Vul?dV<o. 
é 


Then there exists a constant vector Uy such that u(a)>Uy as > oo, and u(x) 
admits, for any surface X enclosing the complement of &, the representation (4), (4’). 

The critical new step in the proof of Theorem 2 is the application of a lemma 
of PAYNE & WEINBERGER [4], which provides an estimate for integral means at 
infinity of functions which possess bounded Dirichlet integrals. 

One might expect that the representation (4) could be used to obtain precise 
asymptotic estimates on the behavior of the solution at infinity. The author is 
as yet unable to do this, the chief technical difficulty being that the components 
of the fundamental solution tensor (9) are not square integrable over the entire 
space. Our best result is the following: 


Theorem 3. Let u(x) denote a solution of (1) exterior to X' and having finite 
Dirichlet integral, {|Vu|?dV <co. Let uy be the constant vector whose existence is 
assured by Theorem 2, and suppose that |u(a)—Uy)|<Cla|~%, «>4, as r>ov, 
for some constant C. Then there is a constant C such that |u(ax)—u)|<Cla|2 
AS L—> Co, 

The same method yields also an improved estimate for |w(a) —uy| outside 
a limited “‘wake”’ region, but we shall not pursue this matter here. 

From a conceptual point of view, it would probably be best to prove first 
Theorem 2. Theorem 1 follows easily from this result and from the work of 
Leray. However, the proof of Theorem 2 requires ideas which are unnecessary 
for the existence theorem, and we shall therefore provide first, in Sections 2, 3, 
and 4, the proof of Theorem 1. We prove Theorem 2 in Section 5, making reference 
to the developments of the earlier sections. Section 6 is devoted to a proof of 
the final theorem, which follows easily from the integral representation (4). 

Throughout this paper the symbols C and K will be used to denote positive 
constants, the values of which may change even within a given context. Thus, 
from the inequality ySC(R?+-1)? we may conclude ySCR for R=1. Similarly, 
the symbol ¢(R) will denote a quantity which tends to zero as Ro. 

Finally we remark that the methods of this paper do not seem to suffice to 
prove the corresponding results for two dimensional solutions of (1). In view 
of the known singular behavior at infinity of the solutions of the Stokes lineariza- 
tion of (1) (cf. [5]), this does not seem surprising. Although this singularity does 
not appear in solutions of the Oseen linearized equations, there may be cases in 
which the non-linearity in (4) prevents the solution from tending to its limit 
quickly enough to permit the Oseen linearization to be a valid approximation 
at infinity. 

2. The Integral Representation 

Consider a divergence-free vector field uw (a) = (u,, “2, 43) and a corresponding 

scalar f(a). We define the second order tensor Tw by the relation 


(Taj =— 6 by +(5, 


4 a 
6;;= , : y 
0, t=—7. 


Ou; a ath 
1 


(5) 


26% 


384 ROBERT FINN: 


Let us also define Lu=Au—Dp. 


We then find easily the identities, valid for any divergence-free fields u and v 
defined in a region V with sufficiently smooth boundary S, 


a) fv-LudV=—Hf (defv-defu)dV+fuo-TudS, 
v v s 


b) f(v-Lu—w-Lv)dV=f(e-Tu—u-To)ds. 


V Ss 


These identities are to be understood in the sense that 


OU; Ou; \ [ Ov; as 
: = -~+ — — OC: Tut aol wy se 
(def w - def u) eS alee Pens)? (LU); 50; 


(6) 


nN = (1, %, 23) = unit exterior directed normal on S, with summation extended 
over repeated indices. Let now X(a#, y) =(X;;), P(«, y) =(P) denote a 3x3 
matrix and 3-vector, respectively, which become singular at x —=y in such a 


way that ax, OX. 
6,=lim [ {R6,— (G4 +) hn as, 
Zy 


7 


where 2, denotes the surface of a sphere of radius 7 about # as center. If, in 
addition, the row vectors of X are divergence free, we obtain from (Ob) the 
representation 


(7) u(x) = — f K-Lu—w-LX)dv,—f (u-TX—X-Tuds, 
Vv s 
valid at interior points of V for any divergence-free vector field u(a) in V. It 


will be convenient for us to choose as quantities (XK, P) a fundamental solution 


tensor for the adjoint system corresponding to the differential equations of OSEEN, 
Au—u,-Vu—Vp=0, 
div ==... 


(8) 


Such a tensor has been given explicitly by OsEEN. It is obtained from the 
relations (cf. [3], p. 34), 
(cf [3] [De 3: ) oo a2@ 


Ox; Ox; 
a) 
f= — [AS —4,-(V®)], 
(9) io 
@ a ; da, 
0 
_ | %| = Moi AY = 2) = 
ey, Ss y+ a [ ut] ) r=|y—2|. 


We have here neglected an inessential constant factor. 


This tensor has the property, in addition to those mentioned above, that for 
every fixed 7, the vector X, ; and corresponding scalar FP; considered as functions 
of y, define a solution of the adjoint system 


Au+uy:-Vu— Vp=0, 


divu =0 


(8°) 


whenever a + y. 
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We now rewrite the equations (1) in the form 
A(u — &) — u)-Viu — u&) —Vp=(u—wW)-V(u— uw), 


div (u — Uy) = 0. 


Throughout this paper we shall fix attention on the vector u — uy, which, changing 
notation, we shall write simply as u. Further, we make the normalization?, 
inessential to the mathematical development, that uy=(1, 0, 0). The equations 
(1) transform to 


ou 
Au Vp=—u-Vu, 
(1’) ot 
divu=—0O, 


and we find the integral representation, valid for every solution of (1’), 


Ox, 


(10) w(x) = AOC KY at Pua Pu KGS! 
10) U(x Le + J u-Vu +f u—u )d 


We choose for S the annular region &y bounded by the given boundary » and 
by a sphere 2’y of radius N sufficiently large that X lies within 2. We consider 
in &y a solution uy (a) of (1’) such that w,(a#) assumes the prescribed values 
u* on 2’, and such that uy (a) =0 on Ly. From the work of LeRay [2], we are 
assured of the existence of uy (a), and we know that the Dirichlet integral 
Dy [Uy|<4 for all N. Since the integrand in the first term on the right of 
(10) is a divergence, we find 
Un (x) = —f xX. (Uy > Vuy) dV + 

En 
ay) Xe Tues pesky d Sf XT, dS, 

D3 “vy 


where the symbol 7, denotes the x, component of a unit normal to 2. 
We note, for later reference, some properties of X and of P. Recalling the 
normalization U)= (1, 0, 0), we find from (9), 


1) Rede in a neighborhood of =O, and also in the paraboloidal 
Je 


? 


“wake” region defined by 7+ (y,— x,) Sconst. 
ii) |X;;|< a for r+(y,—%,)2Clogr, r1>o. 


(12) iii) for all derivatives of order m of the components of X, we have always 


A pean Og ak 
iv) pia y—>0 and r>o. 
Y 


v) Fora sphere of radius 7 about #, we have di Xi (we, y) dSy<Crt, r—> oo. 


aay, 


2 This normalization can be achieved only in the case uy + 0. If uy= 0, the behavior 
of the fundamental solution tensor at infinity is quite different from that of the 
tensor we consider. Our proof, however, remains valid after modification of minor 
details, and provides a new demonstration, in some respects simpler than that of 
Leray, that his solution w(a) tends to zero at infinity. The rapid decay of the com- 
ponents of X outside the limited ‘‘wake”’ region in the case uy + 0 is an essential feature 


of the present paper. 
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From (3) one concludes, by methods analogous to the usual methods of poten- 
tial theory (cf. [2]), that for a suitable sequence of values N + 00, Uy (x), together 
with its derivatives of first order, converges uniformly in any fixed compact 
subregion of €UX. As Noo along this sequence, the integral over 2 in (11) 
tends to a solution of (8) and by (42), this solution +0 as 7 oo. It will be shown 
below that the integral over &y in (11) is absolutely convergent, uniformly in NV. 
For fixed x, it tends to a limit as N-> w along a suitable subsequence, but this 
limit may not be the integral of the corresponding expression involving the limit 
vector. However, the integral of the limit vector exists and tends to zero as 
a—> co, and we shall be able to use this fact. To do so, it will be helpful to write 
(11) in a more tractable form. Observe that,the fundamental solution tensor is 
not unique, and remains a fundamental solution tensor under the addition of 
an arbitrary solution tensor. The integral representation (11) thus remains valid 
if we replace the tensor X (a, y) by the tensor 


(13) X (x,y) =X (x,y) — X(0,y) 

provided the origin of coordinates is chosen interior to one of the finite regions 
bounded by X. The tensor X has the advantage that certain types of singular 
behavior which may occur in integrals involving X (a, y) at large distances from 
» are cancelled by the corresponding integrals involving X(0, y). 

A further difficulty appears in the study of the integral over 2’). Our proce- 
dure here will be to obtain a preliminary estimate using known information 
about the other terms in (11), and then to deform a modification of this integral 
into a suitable volume integral in such a way as to permit a satisfactory a priori 
estimate as N-> oo, 

3. Preliminary Lemmas 


In order not to interrupt the continuity of the proof, we establish here some 
elementary preliminary results which we shall apply later. 
Lemma 1. Let u(a) be continuously differentiable in the closure of Ey, u(a#)=0 
on X'y, and suppose f|Vu|2?dV <K. Then there is a constant K,, depending only 
6 


N 
on K (and not on x or on N), such that 


The inequality is essentially that of Porncar& [6]. The proof we give is due 
to Leray [2]. Under the assumptions, we may extend the region of definition 
of the vectors u (a) to the entire interiors Sy of X\y, in sucha way that the Dirichlet 
integral remains less than a fixed bound. From the identity 


[ t4eav=— f[ V7-Ve) av+ f{ i3&as 


On 


Sy Sy D>) 
with f=w=uj, g=4 log (a —y)?=1 log v2, we find 
U Om; Xx. 
, i Gy, , 
un [pv=-2 [ A avaant] paar 
N Se Ss: 


from which the result follows directly. 
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Lemma 2. Let v(x) denote a solution of the auxiliary system (equations of 
STOKES) 


(15) Av—Vp=0, divw=o 


in the interior of a sphere Xn of radius N, and suppose that on Ly, |v| <KN-*+4, 
v,|< KN~*—3, where v, denotes an arourary derivative of v in a direction tangential 
to the surface =e Then | v («)|< KN~*** throughout the interior S x of Xv. Interior 
to a concentric sphere of pa N/2,|Vv|+|p|<KN~*"4, provided p is normalized 
to be (say) zero at the origin. 


Proof. The vector 0 (#) =v(N«x) defines a solution of (15) in the unit sphere 
»,. Consider the tensor and vector defined by 


i An) Opa Vi 45) (Vp 44) 
Gi; 8x e ! ye cae | 

fi 1 Vi — XG at 
Levis at ees 


Let Yj denote the tensor whose column vectors, as functions of y, define solu- 
tions of (15) interior to 2, and are equal to the column vectors of GY; on the 
boundary. The existence of Yj; is assured by the work of Opgvist [Z]. Let F”’ 
denote the corresponding pressures. The ‘‘Green’s tensor’’ G (a, y) determined by 


(16) G,; = G;, ;- Gi, P, =F! — BP 


then defines a fundamental solution tensor in 2 such that each component G,; 
vanishes on the boundary’. We have the estimate of Opgvist, |TG|<Kr 
uniformly in the closure of the sphere. 


The representation (6b) gives 


For a point # =£0, we find 


5 (a) — af. =|= flowy) —v(N ap 794s 


| | 


since every constant vector is again a solution of (15). We now divide 2, into 
a hemisphere o, which contains the point «/|a| and is symmetric with respect 
to this point, and into the remaining hemisphere o,. We find, in a suitable 
system of spherical coordinates, 


Pe) 5 ey) Sen J elPelsine deo pe: ates, 


where we have used the mean value theorem to estimate E (Ny) —v(N ell 


But |TG|< aa, and |sin g|<|a—y]|. Hence, the first term on the right is 
¥y 


bounded by KN? % The second term is clearly bounded by KN¥?~*, and the 


3 This tensor has been determined explicitly by OSEEN ([3], p. 105). The estimate 
|7G| <Kr-* can be based alternatively on the formula of OsEEN. 
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result on follows directly. The estimate on |Vv|+|p| is obtained immediately 
from the representation, since 7@ remains uniformly bounded when a lies 


interior to a sphere of radius $. 


Lemma 3. Let F'(§) denote a vector-valued function which is defined for all 
points § =(&,, &, &) which lie in an open connected set and have rational co- 
ordinates, and suppose that for some integer n>0, all difference quotients A” F (&) 
of order n vanish identically for all differences taken with respect to such povnts. 
Then F(&) ts a polynomial of degree at most n—1, 

(17) Fé) = s Cr eat 

itj+k=n—1 
where the €;;, are vector-valued constants, and the terms not indicated explicitly 
constitute a polynomial of degree n —2 at most. 


Proof. It clearly suffices to prove the result for a cube contained in the region 
of definition and with faces orthogonal to the coordinate axes. Consider in such 
a cube a line segment L: a<é,<), &,, £,= rational const., and choose on this 
segment equally spaced points 2%, %,=a%)+h,...,%,_,;=%+nh, having 
rational coordinates. Denote by P(é) the polynomial of degree » —1 such that 
P(t)=F(@) for t=a%,...,%,_,. On L set G(t)=F—P. Then G(é)=0 for 
t=2X,...,%,-,, and A”G(é) =0 for all difference quotients of order n taken 
with respect to rational values of t on L. If G(t)=+0, then there is a rational 
point x, on L such that G(x,) +0. Choose =a, and note that for this t succes- 
sive difference quotients of order <m—1 and increment h vanish, since G(E) 
then vanishes at all points of evaluation. Using the point #, to form a final 
difference we find that for this choice, A”) G +0, a contradiction. Hence G(#) =0, 
and we conclude that F(§) is a polynomial in &, of degree not more than n —1, 
the coefficients of which may be functions of &,, . Consider now the coefficient 
of &"*. If there were a difference quotient of this coefficient which did not 
vanish, it would again be possible to form a non-vanishing quotient of F(&) of 
order . Hence this coefficient is constant. Proceeding in this way we obtain 
after a finite number of steps the required form for F(é), q.e.d. 


4. Proof of Theorem 1 


For the annular region bounded by XY and by 2'y, we have from (41) and (6) 
the integral representation 


Uy (x) = —fX- (wy: Vuy) dV + 
(18) sie = — 
2H (AO te el Wy Kn) AS [X- TuydS, 


where X is defined by (13). We know from the results of LERAY that f|Vuy|?4V<K 
& 


independent of N, and that the sequence of integers N may be = chosen that 
Wy converges uniformly together with its derivatives of first order on 2’ and in 
every compact subregion of & to a solution of (1’). We restrict attention to this 
sequence of integers and consider the limit function u (x). Clearly f |Vul?dV<K 


for any fixed compact subregion G, hence, f | Vul?dV<K. Choose e+0. Choose 
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N) so large that f |Vul?dV<e?. Choose R>N, and so large that 
€—-ENno 


(i) i <e, and so that when |a|= R, 
(ii) | f(w-TX—X.Tu—X.-u)dS|<e and 
(iii) | f X-(w- Vu) av| <e. 

ENo 


These inequalities can be achieved because of the properties (12) of X and of P. 
Further, it is clear that for all N sufficiently large (i), (ii), and (iii) will be achieved 
also when w(«) is replaced by the approximating vectors Uy. 


We now consider a value of # such that |a|=R, and write (18) in the form 


Uy (x) = f X(0,y) - (Uy: Vuy) dv, — ad tu X,Y) - (Uy: Vuy) dv, — 


EN, 

— f X(a,y)-(uy-Vuy) i, — f X(a,y)- (uy: Vuy) avy + 
@n,—ENo Gn-€6N, 

+ f (wy: TX(0,y) + wy X(0,y) m — X(0,y)- Tuy] dS, — 
Ps 

— f [uy TX, y) + uy X (a, y) m4 — X(w,y) - Tuy] dSy + 


os Me Died Se 


=04 Qe OF Oss O 0), 


where we have chosen N,=R%. With the above choice of Nj), R, N, we have 

|IO—C,|\<Ke, =f X(0,y)-(u- Vu) dV,. This integral is aug conver- 
e 

gent, for it is dominated by |. f X20, y) u(y) dV, uJ! Vul? dv, 


follows directly from Lemma 1, from (12), and fot the choice of N, and N,>N,. 
We have also |(@)|<e by (iii) above. Further, 


“ail the estimate 


\OPsfKuydVy, f |\Vuyl?dysk e 
: En 6Ni—ENo 
for all sufficiently large N, by Lemma 4, (12), and the ner of No. Lo study 
, observe that for any component of X, =|X (x, y) —X(0, y)| =|X'(E y)|R, 
se & denotes a point on the segment joining the origin to the point #, and 
the derivative X’ is taken in the direction of this segment. Thus, for yc & — &y,, 
N,=R?, we have, using (12), 


ee eee 
vy? Rey te 


by the property (i) of R above, where r=|x—y|. Hence, 


|G Peep “Nav, { |Puy|*dvys Kee 
ON En 
for all N>N,. 
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Clearly © tends to a constant, @=>C,, as NS cor (By? the choice of R, 
|®|<e. Therefore, 
(20) Un (x) =C+e€y(X)+@ 
where €y (a) is such that Jim |€,(«)| tends to zero as R-+>o. 
—>0o 


Now let N->co for fixed a. The left side of (20) tends to a limit, and all 
quantities on the right except (7) have been shown to remain bounded. Hence, 
also @ is bounded. We may clearly choose a subsequence of values N-> co 
such that for any point # with rational coordinates, all quantities in (20) tend 
to finite limits as N-> oo. For this sequence {N} we obtain a function 
im f X(a,y)-Tuy(y)dSy 


Fit) 
( ) N—-oo Sy 
defined and finite for all rational a of &. Corresponding to arbitrary rational 
differences, we may form the difference quotients 
AMF (a) = lim f AMX (ax, y) - Tuy (y) dS,. 

Soa 
Since the row and column vectors of X are divergence free, and since the elements 
of X are majorized by C/r in the neighborhood of the singular point, we have 


AMX .dS,=0. 
oN 
Hence, by the work of Opgvist [/] we are assured of the existence interior to 
Ly of a solution v@ of the Stokes equations (15) such that v”=A™ X on Ly. 
Since Uy (x) =O on 2'y, the identity (6b) gives 


-s 


UN AV 
Ox, iia 


{Bre TuydS = [ vf. (uy-Vuty) dV + f oo 
én ; 


2H EN 


“ap f (Uy: Tv ae ve) -Tuy|dS. 


From Lemma 2 follows directly that the integral over XY tends to zero as N+ 
for any n=O. But 

| fee Pay) av] SK f (0)? ad av, 

én én 

| few a |S K if (o)2aV, 
én én 

by the Schwarz inequality. By Lemma 2 and the property (12111) of X, if ~ is 
sufficiently large, the right hand integrals tend to zero as No. (For our 
case it is sufficient to choose n=2.) Hence, for all larger than this fixed 
integer, we have A“ F(a) =0 at all rational points x of &. Applying Lemma 3 
we see that F(a) is at the rational points a polynomial P(e) m4 5 Ha Bas But 
Uy (x) >U (a) at all points exterior to XY. Thus, from (20) and from the continuity 
of u(x) we obtain that for any a, the limit function u (a) differs from P (ax) +C 
by at most a quantity which becomes arbitrarily small as a—> oo. On the other 
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hand, we know that 
[ 3 aV <0 
x Y 
é 
(Lemma 1), and this can only be if P(#)+C=0. Thus, lim u(«) =0, 7.c., the 
velocity vector tends to its prescribed limit at infinity. “~~ 


To prove (4), choose 
hee (tora |= 


and observe that we can now write (18) in the form 

Uy (x) + f X- (uy-Vuy) dV — ak X- Tuy —uy-TX—uy-Xn,) dS} < ex(R,), 
CR, 

where im _€y (By) +0 as R,-> oo. The desired representation follows by letting 


N-> oo, then R,— oo, and reverting to original notation. To obtain (4’) we begin 
with the representation (cf. [3)) 


py(@) = —Jp-Uy Vuy)dV+ f (p-Tuy—uy: TY —uy- pn} d5, 


EN 2Z+2y 


where we have set pb = V(—), vy =|a—y| and introduced a corresponding pres- 


sue. g = — — (—}. To estimate the volume integral, we again choose Ny so 
Al 
that [|Vu|?dV<e, and then choose R so large that for|w|=R,| fp(u-Vu)dV|<e. 
E€—ENo ENo 
In the region &y— éy, Ny=2R, we have 
fy (uy Vy) avi ss < rel 4 flay \2dV 


which ->0 as Roo uniformly in N (N eae large). We cannot estimate 
the remaining integral over &y — éy, in this way because of the large singularity 
of p at y=0. But we have already proved that u(a#)—>0 as woo, and from 
this follows easily that Vu(a) +0 as a co (cf. [2]). Thus, for sufficiently large 
| x], we shall have | wy (a)| and | wy («)| as small as we wish provided only that NV 
is large enough (depending on |a|). Hence, by choosing R suitably large and 
considering a unit sphere S, centered at a, = R, we can have 


| f(y: Vay) av| <é 
Sy 
for all sufficiently large N. Finally, 


if (wy: uy) av les I IZ | \[Vul2aV 


€—EN> 


ENn,—ENo—S1 
< re 


by the choice of No, for all large enough N. It follows that the volume integral 
becomes arbitrarily small for large ||, uniformly in N as Nov. 


The corresponding estimate for the integrals over 2’ is obvious. 
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All integrals over X'y tend to zero as Noo, & fixed, except possibly 
[w-TuydS. Hence, this integral remains bounded as N-+oo, and we can 


=v 


choose a subsequence of the {N} so that 


F(x)=lim fp-TuydS 
N— oo hy 
is defined and finite at all rational points. Taking a single difference quotient 
A F(a) and repeating the previous procedure of replacing the difference quotient 
AW by the solution v@) of (15) equal to J on Xy, we obtain an estimate* 
on v sufficiently good to ensure that AY F(x) =0, 7.e., M(x) =const.= Po. 
Inserting this constant in the above representation for p(a#) gives (4’), and com- 
pletes the proof of Theorem 1. 


5. Proof of Theorem 2. Alternative Proof of Theorem 1 


We base the proof of Theorem 2 on a lemma of PAYNE & WEINBERGER [4], 
concerning integral means of functions with finite Dirichlet integrals in a neigh- 
borhood of infinity. We shall need, actually, a somewhat more general result, 
for which we provide a short proof. 


Lemma 1’. Let u(x) denote a vector-valued function with continuous first 
derivatives, defined in the exterior Ep of a sphere Xp, of radius Ry, and suppose 
that the Dirichlet integral Ap [u| extended over Ep, 1s finite. Then 


i) The integral means U of u(a) over concentric spheres of radius R tend to a 
hmit Uy as R-> ov, 1.¢., of 


u(R) =), [ u(Ro) dQ, 


Da 


be’ 4 


where w, dQ denote a point and surface element on the unit sphere, respectively, 
then there exists a Uy such that lim u(R) =U. 


li) For any R2E],, 


1 
x | Maw) Ss 4, [u] =| \Vulav. 
=R Ep 
il) If 7,4 denotes the distance between two points x and y of Ep,, then 


2 
Vn y 


f en —wel d ee IK << CS 
So 


independent of x, where K depends only on Ap [u). 


The proofs of i) and ii) are given in [4] and will not be repeated here. To 
prove ii1) we note that w(a) can be extended to the entire space as a continuous 


. 4 The magnitude of the vector field (a, y) decreases sufficiently rapidly at 
infinity that a formal manipulation would imply F(a) =0. However, this vector field 
has the character of a potential flow due to a source at the point y= «x, and creates 
a net flux across the boundary of any region enclosing the source. Therefore, there 
cannot exist a non-singular, divergence-free vector field wy (a) equal to w(x) on the 


boundary. This behavior must be expected, for the pressure is defined by a solution 
of (1) only up to an additive constant. 
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function with bounded Dirichlet integral, depending only on A r, U4]. Consider 
now a sphere Sz of radius RK and center at #. In this sphere the relation (14) 
becomes, noting that w(x) no longer need vanish on the surface Xp, 


y2 
SR Sr 


[Seay =—2 f (us — me) 8 ave t fu —u)eas, 
ny = } : 
=R 


and the result follows directly from Schwarz’ inequality and Lemma 1’, ii). 


To prove Theorem 2, we consider, as before, the annular region &y bounded 
by a sphere Sy and 2. We assume again® that u)=(1, 0, 0), and we write the 
equations in terms of the variable u—wu,, which, changing notation as before, 
we write simply as u(x). We note that since we are dealing with a single vector 
field rather than a sequence, we shall not be bothered by questions of uniform 
convergence, and it is therefore no longer necessary to introduce the tensor 
X(0,y). We obtain the integral representation 


u(x) =f X-(u-VuydV+f (X-Tu—u-TX—u-Xn)dS+ 
En = 


ios) (X-Tu—u-TX—u-Xn,)dS. 
SI 
Certain new integrals over 2’y appear here, which did not occur in (48), 
since u(a) need no longer vanish on this surface. 


The integral over 2’ obviously tends to zero as x-> co. To study the integral 
over 6y, we observe that for any e¢>0, we may choose WN, so large that 


J |\Vul?¢V <e?, and then |x| so large that | fX-(u-Vu)dV| <e. For this a, 
€—-EN> EN 


| f X-(a-Pv) av) s fx wav f \Vul2av 
|\Ev—Eno é €—-EN, 
and hence, by Lemma 1’ and the properties (12) of X, this integral tends to zero 
aS &—> oo. 
We are left with three integrals over 2'y. Two of them can be estimated 
directly, for 


fuTXaslsk fas f Neds, 
Ey se 


2y 


[wXmas = Nel assf Kids. 
iy Sa Sih; 


Each of these terms tends to zero as N->oo, because of (12), Lemma 1’, and 
the fact that the Dirichlet integral of w(a) exterior to 2’y becomes arbitrarily 
small if N is large enough. 


We estimate the final term exactly as we have estimated the corresponding 
term in the proof of Theorem 1, a simplification appearing in the fact that the 


5 The exceptional case u)=0 requires only minor technical changes in the proof; 
cf. footnote ?. 
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limit of this integral is now seen to yield a function 


F(x)=lim f[X-TudS 
N-> oo Sy 
which is defined and continuous at all points, not just the rational points. Again 
forming difference quotients and replacing A™X interior to 2’, by the solution 


v” of (15) which equals A”) X on 2'y, we find 
n —_ (n) . Ou y ; ry 
fa X-TudS=| vf + u Vujav a 
=N On 


+ f [w To — of). TuldS + fu- Topas. 
Pa Zy 


The estimates proceed as in the study of (19), @, except for the appearance of 
a new boundary integral fu-Tv) dS. It follows from Lemma 1’ that this 


integral tends to zero as N ex oo, provided we know that for m sufficiently large 
| Tv|<KN-! on Ly. The proof of this fact can be obtained in a manner 
entirely analogous to the proof of the corresponding estimate of potential theory, 
using the Green’s tensor (16) in place of the Green’s function of potential theory, 
and we omit the details. 


We now apply Theorem 2 to obtain a second proof of Theorem 1. We observe, 
using the original notation, that the approximating solutions uy (x) of LERAy tend 


to a solution w(x) in & for which [ ee dV<co, and for which 
om 


J \Vul? dV <oo. Hence, by Theorem 2, there exists a vector uw, such that u(a)>uU, 
é 


as #—> oo. It follows immediately that u,=Up, g.e.d. 


6. Proof of Theorem 3 


Once more, we replace the symbol w(a#)— wy) by the symbol u(x). We are 
given |u(a)|<Cr*, «>. We have, by Theorem 2, the representation 


u(x) =f (X-Tu—u-TX—u-Xn) dS —f[X.(u-Vu) av. 
5 e 

Since divu=0, u-Vu=div wu, and we may integrate the volume integral by 
parts to obtain 


y 


u(x) =f (X-Tu—u-TX—u-.Xn,—X-uu)dS +f (uu) -VXav. 
é 


We have here estimated a boundary integral as in the proof of Theorem 2, and 
passed to a limit. 


Using (12), we see that the integral over XY has the order r+ at infinity. To 
estimate the volume integral, we choose a point # with distance 2R from the 


origin (K sufficiently large), and divide the region of integration into three parts, 
as follows: 


I: the annular region &% bounded by XY and by a sphere of radius R about 
the origin, 
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II: the region Gp consisting of those points interior to a sphere of radius 3R 
about # as center, which do not lie in &p, 


III: the region Ez of points whose distances from # exceed 3 R. 
We have 

f (wu)-VXav| <fuwtav f |7X|2av 

ER ER ER 


<C(14+R)-4*+8 f |WX|2dV 
AR 


by assumption, where Ap denotes the annular region bounded by spheres of 
radii R and 3R about «#. The integral over Ag can be estimated by a tedious 
but formal computation, using the definition (9) of X (it is necessary here to 
use a more careful estimate of X than is indicated in (12)). We obtain 


f [VXPdvVs CIR. 
A 
30 that ; 


if (wu) VX av| = C(R1+ R-22+3) 
ER 


for a suitable constant C. 
In Gp, we have everywhere | ee] ZAG IR, Soy lien 


J (wu) VXav| <CR**{|7X|av. 
The integral on the right can again be estimated by formal computation, which 


yields 


{ (uu) -VXdV|<CR*+4, 


GR 


Finally, 
f (wu) -VXdV 


ER 


<C fr-**dr f |\VX| ads 
3R Hie 
= C R7-24+3 


by essentially the same calculation. Collecting these estimates, we see that 
| w(x)|SC(R-***44 R4) at infinity whenever «>. Set a=3+2e, €>0. We 
have proved that the inequality |w(a)|<CR~*~* implies always the inequality 
|u(a)|<C(R-#~*-+R4). A finite number of iterations then yields the desired 
estimate |u(x)|<CR4, ¢.e.d. 


Remarks. The estimate of Theorem 3 is exactly the order of magnitude ex- 
hibited by the solutions of the linearized equations (8) in the “wake” region 
7—x,<C. For paths tending to infinity outside this region a better estimate 
can be obtained, for contributions of the order 7+ appear only from the integral 
over the region I and the integral over 2, and these contributions can be shown 
to be of smaller order when the point # lies outside the wake. The nature of the 
representation suggests also the possibility of examining the qualitative behavior 
of a one parameter family of solutions of (1) depending (say) on the viscosity, 
as the viscosity tends to a zero or infinite limit. The author hopes to return to 
these matters at a later time. 
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A. Study of the Interior Ballistic Equations 


RAIMOND A. STRUBLE 


Communicated by C. TRUESDELL 


Summary 


A system of non-linear differential equations, typical of the interior ballistics 
equations for orthodox guns, is studied. Characteristics of the mathematical 
pressure curves are obtained for the entire spectrums of the ballistic parameter, 
burning exponent, and geometric parameter (characteristic of the form function). 
The mathematical solutions are found to distinguish sharply between regressive 
and progressive burning geometries, and a constant pressure solution is obtained 
(a piezometric efficiency of one) for what is defined as the ideally progressive 
geometry. The approximate validity of isothermal solutions, established by 
CLEMMow for a linear burning law, is extended to the general case. 


1. Introduction 


As a technology, the interior ballistics of guns is in an advanced state. Indeed, 
a recent paper [/| suggests that the main problems have been almost completely 
solved. Still, throughout the many contributions to theory, little attention is 
given to mathematical aspects of the problems, and, as a result, the full impact 
of the equations is not always felt. Even in the most orthodox cases the distinc- 
tion between practical ballistics and research ballistics is obscured by a common 
empirical flavor, and useful aspects of the mathematics appear to be lost in a 
zeal for correlation with experimental data. Thus, despite a long and continuous 
history, many fundamental questions concerning the gun equations have remained 
unanswered, and so perhaps a few purely mathematical considerations are ap- 
propriate notwithstanding the rather advanced state of the technology. 

In the present paper we consider the interior ballistics problem as characterized 
by a system of non-linear differential equations. These equations are not unlike 
those appearing in numerous other investigations (cf. [1], [2], [3], [4], [5]). Here, 
however, we are concerned with the implications of the equations themselves, 
so that questions concerning the degree to which these equations or their solutions 
represent the physical situation are not considered. On the other hand, much 
of the mathematical study is motivated by practical aspects of the ballistics 
problem, and some of the results may prove useful either in design problems or 
as aids to an understanding of the fundamental ballistic processes. A list of 
nomenclature is provided as a ready link to practical applications of the results. 
By-and-large the results are new in the sense that they are here derived for the 
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first time from the equations, although some may well reflect established tech- 
nology. Several results, peculiar to propellants with progressive burning charac- 
teristics, are new in any sense and should prove both of practical and of theoretical 
interest. 

Specifically the equations considered are introduced in § 2 and their solution 
interpreted geometrically in § 3. Characteristics of the pressure curves are derived 
in §§4 and 5, and parameter effects are studied §§6, 7, and 8. In Oe 
discuss the ideally progressive case which admits a constant pressure solution, 
and in §10 we consider the significance of isothermal solutions. 


2. The fundamental equations 
The interior ballistics problem to be considered is in reference to the following 
equations *: 
Momentum Equation 


dian 
M— SA, (1) 
Burning Law 
Gf De 
gps Ope (2) 
Equation of State 
NF = PAx (3) 
Ty 
Energy Equation 
NF if 1 dx \2 
faa! ( n= aa) 4) 
Form Function 
NeaG hai ys (5) 


The principal nomenclature used is as follows: 


A =Gun bore area 
B =Propellant burning constant 
C = Propellant charge 
F =Propellant force 
/ =Fraction of propellant web remaining 
M = Mass of shot 
N =Fraction of charge burned 
P =Instantaneous chamber pressure 
Fi, = Average pressure, defined by equation (24) 
P,, = Mathematical maximum of P 
J pea= ee = Closed chamber pressure 
Ai) 
B°MFC je 


W2 A2 


p —(B2MEC\s 

2. ee 
| W Ae 

P, SNAG) ca 


1 
Ballistic parameters 


| 


W?A3 x, 


* See [5, Chapter V], for example. 


| 
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aa 1 
— 3-20 
S =Surface, bounding trajectories, defined by equation (19) 
T =Instantaneous gas temperature 
T) = Adiabatic combustion temperature 
Pe lime 
v =Instantaneous velocity (speed) of shot 
W = Propellant web size 
x ==Shot travel 
%q = Initial position of shot, defined as the initial chamber volume divided by 
the bore area A 


a == Burning exponent, 0<aX1 


B vay? J 


2—a 
y = Ratio of specific heats of propellant gases 


é, — Instantaneous ballistic efficiency, defined by equation (75) 


Te : : ns 
é = —* = Instantaneous piezometric efficiency 
enpeoi 


@ =Propellant geometric parameter, —1<05 +1, 
—1<%8<0 corresponds to progressive burning while 
0<#S +1 corresponds to regressive burning 

we =(—f" 

My, = Value of w at maximum pressure 


y= | ae v = Dimensionless shot velocity 


Ce = Dimensionless shot travel 
0 


Ct=1 =} 


3. Geometrical formulation of the ballistic problem 


Denoting dx/dt and 1—f by v and o respectively, the fundamental equations 
(1), (2), (3), (4), and (5) may be written as the autonomous system 


dx _y 
cepa 
du A 
di M7? ) 
Diem ue pe 
Cie Bee? 
wherein 

| J 1 2 

pa lFatt9 bo)—+s(y 1) Mv? (7) 

A 4% 


A solution of system (6) is interpreted as a space curve in (¥%, v, o)-space with 
time ¢ a parameter along the curve (see Fig. 1). The (%, v, a)-space is called the 
phase-space of the system (6) and a solution is called a trajectory. The parameter ¢ 
may be eliminated by dividing the first and second equations of (6) by the third. 


27a 
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Hence the trajectories are also solutions of the system, 


ax W v 

deo BPP (8) 
dv _ A W pi-« 

Gas Wie IB: ; 


We are interested in those solutions corresponding to a zero shot start pressure, 
i.e., P =O initially. The appropriate initial conditions are therefore 


i= Kg 0; 10 o=0. (9) 


In the phase-space, these initial conditions are represented by a half-line, the 
positive x-axis. Since at each point of the half-line the mght hand members in 
(6) vanish simultaneously, the initial 
conditions (9) correspond to critical 
points* of the system. Each of the 
latter constitutes a trivial solution 
(so called point trajectory) of (6) 
which is, of course, of no practical 
interest. However, equations (6) 
admit other non-trivial solutions 
emanating from the critical points. 
The latter therefore are not deter- 
mined by their initial values alone. 
Fig. 1. Phase space trajectory This lack of uniqueness is a result of 

the fractional nature of the ex- 

ponent «. Indeed if «=1, the system (6) is analytic at each point of the half- 
line (9), and only the trivial solution is admitted: This is an interesting point 
since a linear burning law, «=1, is by far the most popular mathematical 
model, and the corresponding solutions are widely applied in practice. In § 6 we 
shall see that the latter violate the initial conditions with reference to the time 
variable although the geometrical conditions (9) in the phase-space are realized. 


For a<1, the second of equations (8) requires that 


U 
an (10) 


as the initial point x = x), v =0, o =O is approached along a trajectory. Thus (7) 
implies that 


— 2 

CF(1+0—60)—7— m+ 

ba Cs ) 2 a CRUsa) (11) 
a Ax A Ky 


so that along a non-trivial solution, the pressure P vanishes as o at the initial 
Domed Ax: 


From the second of equations (8) we obtain 


ree B(Eieie 


, 


2-a M W\o 
* See [6, Chapter II], for example. 


** Generally —1<#, but for later considerations in § 9, we note here that if 


@—=—1, then P vanishes as o? and Ete en 


o A Xo : 
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and in view of (11), consider a change of independent variable to 


so that 


dv oo A apenas 1—o 
du 2—a M B 


A % gy 
at the initial point. Using (12), equations (8) become 


u 


x 
dem al Ay ke jee 
Cte oO 9b ye \ =P ‘ 


4 (13) 
gee Get S nenwafoap TK F 
dn 2-«a M B ne ; 
pwe-% 
At this point it is convenient to introduce the dimensionless variables 
Y = || Ley, 
CE 
(14) 
Ses 
eae: 
For conciseness, we also introduce the symbols 
me CF 
c 8 liven ) 
Mis (15) 
P =(Fpa)* 
ae Ae 


Each of the latter has the dimensions of a pressure. The first is, in fact, the 
closed chamber pressure, 7.e., the maximum pressure which would be produced 
upon burning the complete charge in the initial chamber volume. We shall see 
in §5 that the quantity P, is closely associated with the maximum pressure of 
the actual ballistic process. If (14) and (415) are introduced into (13), there is 
obtained the system 


dé =p(2 1a /P yB\e y 
0B. ‘i 
dv =6(E) "(ga)" 
dM B) Fy uh 
wherein B = and 
Deh : 
V— Ay 
see ae (17) 
Rye E 
The initial conditions appropriate to equations (16) become 
é=1, 
) == (0), 
v (18) 
Toma 
TP 
ie tea 
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for 4=0. The latter two follow from (10) and (41) using (12). Owing to the 
condition P/P.u°=1-+0-50, the initial point of the trajectory in (6, v, “)-space 
no longer appears as a critical point. Thus (16) exhibit only the non-trivial solu- 
tions. The system remains non-analytic at the initial point, so that existence 
and uniqueness questions remain. However, the latter are settled by standard 
results from differential equations theory*: There exists a unique solution of 
(16) satisfying the initial conditions (18). Such a solution is the non-trivial solu- 
tion of the original boundary value problem (6), (9). The formulation (16) is 
therefore fundamental to the interior ballistics problem. 


4. General properties of the trajectories 


If a solution of (16) is interpreted as a space curve or trajectory in (6, y, )- 
space, then the trajectory is limited to that portion of the first octant for which 


&=1. In addition, since the solution exists only for oi = 0, equation (17) implies 
c 


that the trajectory must also lie 
within the surface S defined by 


f (440) Pou y—*9?=0. (19) 


For #>0, S is an elliptic-shaped 
cylinder, bounded in the uw and » 
directions (Fig. 2), while for #S0, 
S is a hyperbolic-shaped cylinder, 
Fig 2) Suriace, 5 for 0s open in the mw and »y directions 
(Fig. 3). In addition, a physically 
significant trajectory is limited to the halfspace wX1, since the charge C is 
completely burnt at «=1. However, as this limitation constitutes an artificial 
discontinuity in the system (16) which bears no mathematical significance, 
we shall not impose it 
generally **. 

Figs.2 and 3 suggest 
that the complete mathe- 
matical trajectories for >0 
may differ considerably 
from those for ?<0. That 
this is indeed the case we 
shall see in §5. However, 
some observations of a gen- 
eral nature may be made 
here. All trajectories ema- 

Fig. 3. Surface, S, for 9<o0 nate from the same point, 

&=1, p=7=0; with a 

tangent in the €=1 plane of finite slope. This is in contrast to the vertical 

(v-oriented) tangent at the initial points of the trajectories in the phase-space. 
* See [7, Chapter I], for example. 


Sok 
Bee 5 discussion : the complete space curve solution is both instructive and, as 
we shall see in §§ 7, 8, and 10, of practical concern. 


{ 
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Since a = 0 throughout, a trajectory continuously rises in the y-direction and 


twists toward § increasing. Such a trajectory will thus terminate only as it 
contacts the surface S, at which point the pressure P will vanish. For #>0, 
it is possible for a trajectory to avoid the surface S only if it progresses up (& in- 
creasing) the elliptic-shaped cylinder. However, in such a case since > oo and 
both w and y approach finite limits, (17) implies that P will vanish. Thus in any 
case for #>0, P must ultimately return to zero for some finite value of w and 
the mathematical solution terminates. We note in passing that P necessarily 
attains a mathematical maximum. 


In contrast, we shall show in §5 that for #<0, the pressure cannot return to 
zero so that the mathematical solution may be continued for all « . In addition, 
P will approach zero as w—> oo only for #=0. Hence, a mathematical maximum 
for P when #<0O appears uncertain. It is well known, for example, that some 
solutions for the linear burning rate, « =1, exhibit continuously increasing pres- 
sure characteristics. In §§5 and 6 it will be shown that « = 1 is, however, exceptional 
in this respect. 


5. Pressure characteristics 
We rewrite (17) as 


= (14-8) pF — bur? Poh yp (20) 
and compute the y-derivative of (20). Using the second of equations (16), we obtain 


d phe ae ae , 2 wee fe i Ma) 
a (Es) = (4 0 26u°) (3) (y —1). 
On the other hand we have 


we) aletee 


dé oP : 
7 ase 21) 
Hence, 
eg ee OV vale Ae er a 7 5 
> dy (2) - Ge | yb ee B) yb | (22) 


Considering once again the second of equations (16) expressed in the form 


Nilo. 
t= (Bayt 5 
1 
we may write 
ub 
ee ee 8 
v= f(z) au 


Hence the quantity »/w® is the average value of (P/P)*~* as a function of we. 
It is convenient here to introduce a function P, which depicts an average value 
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of P itself. Accordingly we define 


a-tal-[e fare " 
0 
Thus (22) may be written as 
fate) Ce) 2 Te ae ha ana 


We now observe that P cannot increase without bound, for then the third 
term in the brackets of (25) would dominate the other two and d/dy (P/F) would 
ultimately become negative. Thus the pressure, P, must exhibit a mathematical 
maximum. This was shown in § 4 for # >0 but now applies generally. For #<0, 
(25) ostensibly permits P to increase monotonically to its maximum (finite) 
value as > co. However, the y-derivative of (25) reveals that in such a case 
d?(P/P,)/dy? would ultimately become and remain positive, which in turn would 
lead to an unbounded P/P,. Hence, in every case a maximum*® of P, which 
we denote by P,,, occurs for finite w. The value of the latter we denote by uw,,. 


Since (dP/P)/dyv vanishes for P =P, and certainly P.<P,, (25) requires that 


a (26) 


and a fortiori 
20 Eine 
(z) 


1 


For # >0, the latter is trivial but for ®<0 yields the important inequality 


1 
DO BANVIE EGA more 
a y WwW? A2 | Sn (27) 


An upper bound to F,, is also available. In fact, using (24) we may rewrite 
(17) in the form 


i ah iy ca Set YE a 
iD ole ait 2 (3) 
PR pee aa ; “J 


ia 
S 


Thus P/P,?? decreases monotonically up to maximum P. Rewriting (23) in 


the form 
d eeey.? P a 2—B 
5 rad) ai ae 


we conclude that 


or 


* Strictly speaki . , 
se y speaking, we refer here and throughout to the first mathematical maxi- 
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up to maximum P. Hence from (22) we have 
1+¢6 En, \P-?# 
sp 202 ry (2) ; (28) 


and (26) and (28) together yield 


al 1 
eg mel eden on 


we, 
y 


eae | Mm 


IA 
| 
mn 
| 
bo 
= 
L 


(29) 


iy 


It might be noted that 3!< (3 —2a) *~** Se for 0051. Actually the-bounds 
in (29) are not determined since the value of 3, is in general unknown. However, 


for the special case ? = —1 considered in § 9, (29) becomes 
esas pail fie 
(St 2-56 = 20)2 ea (30) 
iy es Y 


It should also be observed that as the maximum pressure is delayed, 7.e., 3, 
becomes large, (29) approaches the explicit form (30). Another observation of 
interest here is that according to (25) the unknown quantity in (29) is given by 


1+¢0 
= ee see ee (31) 
y F, Te at 
Using (31) in (29) results in 
on 
Pies Page a) Py (32) 
which relates to the piezometric efficiency. For convenience here we define* 
the latter by the ratio ¢,= — so that (32) yields Sa eS 
be (3—2a)i-@ 


In § 4 we observed that for # >0, the pressure P returns to zero for finite ju. 
However, if <0, from (25) we see that P cannot decrease to zero since d(P/P,)/dv 
would then become positive **. 

Thus for #0 it is of some interest to consider the possible hmiting values 
of Pasuw—>oo. It is not difficult to obtain an asymptotic solution of equations (16). 
The latter is a constant pressure solution with 


20) p = 28 SMEG Ver 
ae 1 5 


P=(— 


It is interesting to note that (33) is given by the lower bound to PF, obtained 
previously, (27). 

We infer then that for #0, the pressure curve rises to a maximum which 
exceeds the value in (33) and subsequently decays and approaches the asymptotic 
solution (33). In §9 we shall again consider the latter together with its special 
significance for ideally progressive burning propellants. 


* Piezometric efficiency is generally referred to muzzle values only and depicts 
a displacement-average pressure to maximum pressure ratio. 
xx For @=0, this argument fails with ~— co since P may tend to zero with the 
au) 


remaining term a 
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6. The linear burning law, a=1 


It is well known that ballistic solutions can be obtained in closed form for 
a —=1. Indeed much of ballistic technology is based on these solutions despite 
overwhelming experimental evidence suggesting that the appropriate values of 
« should be less than unity. It is important therefore to understand in what 
sense the «=1 solutions and the «<1 solutions resemble one another. Not so 
well known are many significant differences: Actually, the linear burning case 
is exceptional in many respects as it exhibits numerous phenomena not exhibited 


by any other case. The results in §5 indicate some of these variances. For 
es . —20 B*MFC 
example, from (27), it is to be expected that when the quantity > WPAB 


exceeds unity, the maximum pressure, P,, increases with increasing «. In 


fact, for «1, (27) requires that P,—co. Hence, the »=1 solutions would 


pe 2 
not be expected to exhibit maximum pressures when 2° . ie ota =41.-The 


explicit solutions, which we shall consider presently, verify this conclusion. How- 
ever, the latter is in sharp contrast to the fact that for «<1, every solution 
exhibits a maximum pressure. 


To consider the limit «1, we rewrite equations (16) in the form 


dé _ » Pie (ae » dy «» BI* 7 P \ie 
du P i eer we? dp B P, (za) j G4) 
: = [ESI IE Oso. ; [PB ve 
wherein P,—P} “=| : 
rein f= weaz} 3 independent of «. Since Pb is bounded away 


from zero initially, the solutions depend continuously upon the parameter «. 
In view of the above remarks concerning maximum pressures, we conclude that 
SONU SIMUETEC, a} Pe 

AVA ASS a ya 
pressures for «>4 must occur for increasingly larger values of u as «—>1 (2.65 
Ly o). 


With « =1, the solution of (34) is* 


when the quantity exceeds unity, the maximum 


2 
page 2) Fe, ___ |e Pety—7 
2(1+9) P,— (20 P2+y—1) : 


The corresponding pressure ratio is 


1p 1 


2 
PP 4 Bly Pah nL DE Or aay L hanes 
stander Blom 2 Es Orr) whl core era alee ae aS 
so that 
d (P/f) P 2P, o) — 2 es 
AP /ee) __ 149) —(2.0P2 : ? a(1+9)—(20+y—1) » e@aptpe 
dy a1 + 9) — 20 P+ y) 9] SACHS MOMMA ices ae cae ore? 
* We assume here that 20 P?+y—1+0. If 20 P?+y—1=0, then E=exp.- 4 


(1+8) P, ¢ 


Interior Ballistic Equations 407 


Hence a maximum of P/P. occurs for* 


pte) te 
TRE 9 pasta (38) 
if 0<20P24+y, ie., if Se PE <A. 
The maximum value of P/P, becomes 
20 Paty 
Pe Peay? poe riyenees 
P~ 20PP+y |2(20PR+y) 
and, if 2023 --y—1 >0, ous subsequently returns to zero for y= pri a 
C Dr ia 
while E-+co. However if 208 P?+y—1<0**, returns to zero only as yo. 
a 2 c 
For 20 P?+y<0, 4.6., a28lt > 1, equation (37) implies that ee increases 


monotonically without bound as y>oo and E—ow. 


These results give strong indications that parameter effects on various ballistic 
quantities such as maximum pressure with «<1 may not be adequately reflected 
in the « =1 solutions. Some considerations of these effects appear in §§ 7 and 8. 
Further evidence of the peculiar nature of the «=1 solutions is afforded by 
considerations of the corresponding phase-space trajectories. The latter are 
defined by the equations 


Maes 
Bel (9) 
=a. 


where € and » are given by equations (35). Although x =x), v=0, o=0 is the 
initial point of the trajectories for «<1, it is attained in the limit only for t>— co 
along the limiting trajectory (39) ***. Indeed, equations (34) with « =1 imply that 


(Regd au dy (40) 


Ry ORE Piao ap 
is not integrable at »=0. However dé/y is proportional to the differential of 


time dt. Hence the parameter ¢ cannot be introduced **** along the trajectory (39) 
so as to yield a solution of equations (6) satisfying the initial conditions (9). 


2(1+0) & 
26 Fe+y—-1 
however, exceeds the value in (38), so that the maximum is actually attained. If 
20 P? + y—1<0, then y is not bounded. 
xx This inequality together with the condition 20 PF? +y>0 requires 0<20 P?+ 


Veal: : : 
*xkx The solution (39) might be interpreted as a 3-dimensional separairix [6, 
Chapter V]. 
*xxxx A device such as used by CorNneER [5, Chapter 4] to avoid the singularity 
in (40) at y=O merely truncates the solution so that both » and P fail to vanish 


initially. 


x Note that if 20P?7+y—1>0, then y is bounded by , which, 
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7. The geometric parameter, d 

In this section and in §§ 8 and 9 we consider similarity transformations which 
form the basis for parameter studies*. These transformations are not original 
here, although the present treatment and use is somewhat novel. 

The importance of the geometric parameter # has been noted previously in 
the general discussion of the solutions. In particular, it was observed that the 
two ranges <0 and #>0 exhibit marked differences with reference to the 
trajectories and the pressure curves. Still, the solutions of equations (16) depend 
continuously** upon @ in the range —1<#<1***, so that some similarities 
must also be present. The latter may be exhibited through the similarity trans- 
formations. 

We note that for #=+-—1, equation (17) may be written 


26 Ya ks, y? 


a ee fea err 7° j 
1-29 Pe é ; (44) 
For ? =1, (41) becomes 
pen ae 
Je: ass Dims i 
ae E : (42) 


Equations (41) and (42) suggest the change of variables, say for #>0, to 


- =| 28 ae 
Los ves ) 
(43) 
— 2 V. 
Aa 
For then (41) becomes 
y—1 vi 
40 Bek lee er 
(+P Pu’ E Ce 


which is (42) except for notation. It should be observed that equations (43) 
merely effect scalar changes in variables which are subjected to homogeneous 
initial conditions. Introducing the new variables (43) into equations (16) and 
using (44), we obtain 


Lene (1t+a2 pli E |* » 
a B REO 7 ; 2 rae 
My 4ohi=aa) 1 2— we — vee ae be 

2 pe 

45) 
2 g_ Y—1 WY y1-« ( 
dm _ [+e Rya|? 7 ae 
3— 

diy a eae) 1 &€ ; 


easy 
Ae R6GGLA [4] and CORNER [5] have employed similarity transformations exten- 
Sively as a means of organizing tabular solutions, 7.e., ballistic tables. 


I u 
2 ve 
This statement nust be q alitied of course, with reference to the extent of 


*** § — — 1 is exceptional and is considered in § 9. 
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which are equations (16) in the new variables (£, 4,,»,) with @=1 and P/R 
+92 P 


replaced by ~~ (a=2ay a: Hence any solution with O0<#?<‘1 is given by a cor- 
49 2-20 1 

responding #=1 solution. We note here that the pressure ratio P/P, in (16) 

is replaced by a eeou = in (45). Since € is invariant, it is convenient to regard 


— as the independent variable throughout. Equations (43) are then interpreted 
as functional equations holding identically in § (€ =1), and all properties of the 
solutions for 0<#<1 are reflected in the #=1 solutions*. In particular, the 


pressure ratio is given by 
P_ (149)? /P 
a. ( ie eu 


f AR NE, 


holding identically in €, wherein (P/P,)g_, is the pressure ratio given by (16) 
for ?=1 and for an equivalent ballistic parameter 


The effects of variations of ? on the pressure ratio are therefore two-fold: A simple 
scalar effect as given by (46) and the more subtle effect inferred from # variations 
on the equivalent ballistic parameter, (47). The similarity transformation (43) 
limits # to the range O0<#1. This, of course, is to be expected from the discus- 
sions in §§ 4 and 5, since the solutions for 0<#X1 terminate with finite yu (7.e., 
finite 4), while for —1<8S0, they are defined for all wu. 

For # in the range —1<#<0, similarity transformations are also available. 
By way of illustration we may choose ? = — $ as a standard. The transformation 
to other ? in —1<#<0 is then defined by the equations 


Si eee 

_ 3 48) 
rss 
ne ee 


wherein the triplet (£, v2, W) is the trajectory given by equations (16) for? = — 3, 
and for the equivalent ballistic parameter 


(=) = (1+0)2 < , (40) 


The pressure curve is transformed according to 


P 2(1+8)2 (P 
eo erie ee oe) 


¢ 


1 
z 


* It might be noted, however, that ‘“‘burnt’’ determined by ~=1 now becomes 


z 20 ? 
4 Ue 14+06 . 
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In this instance, however, it should be noted that ‘“‘burnt” is now given by 


1 
p= (=). which varies between 0 and -+-co with #9 in —1<0<0. Hence, 
the complete mathematical trajectory is required if the similarity transformation 
(48) is to be used in general. Again ¢ is invariant, and all properties of the solutions 
for —1<%<0 are reflected in the #=— 4 solutions. 

The similarity transformations (43) and (48) amply demonstrate the con- 
tinuous dependence of the solutions upon # in the two ranges —1<#<0 and 
0<0<1. The two limiting cases #->+0 and #-——1 are not covered, however. 
These are treated (although rather incidentally) in the following two sections. 


8. Ballistic parameters 
We consider a general similarity transformation 


E=€, 
ts=($) & (51) 


wherein @ and 0 are constants. In terms of the new variables, equations (16) 
become 


OS” P, \l-a—B (l—a) a—b 
d bs Wr z) - 
x § + ois 
P.\—84 = —2b+8 1 Fie 
(1+8) wg— ($5) 9 ne— 2 1 (2) “a | we? 
1 1 
Diy P,\1—a—B (2—a) a+b (52) 
d ls a (z) iS 
eae —1 /P,\-25+Ba = 
(1+%) ue— (=) epee (¥) 2 |l-« 
5 3 P. M3 2 PB V3 1 
& i ee 


We should like to choose a and 0 so as to eliminate any effects of variations of 
the ballistic parameter P/P, in (52). This is not possible in general, but if we 
choose a and b so that 

a—B(1—a)a—b=0, 


1—a—B(2—a)a+b=0 


2—20\ 1 — 
(ie. a= (= =) b= aE then equations (52) become* 


a =2/f'| st eevee? *¥s 
a [ls jer (Eee a 2 2 
1+9— Aan 8 ug—2—* e M3 
3 
Sam Staal mere (53) 
ea 3) 3— 20 Ga vam 9 
es ee oe 
d ky &é j 


* Note that rather fortuitiously —2b+ Ba=o. 
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and (51) becomes pee 
M3 = (=) Per Seve 
‘ (54) 
Cl 
Tae [Baltes 
1 
while the pressure ratio is given by 
2-20 (1+0) igo (2) 088) 2B jammy | ye 
Pe P\\gaea) MS ae se ‘ 
plz) oo 


Hence any effects of variations of the ballistic parameter are now reflected either 
in the coefficient of the [—#u§] terms in (53) or the simple scalar effects in (54) 
and (55)*. 

For #=0, the [—@y§] terms are not present, and variations of the ballistic 
parameter can effect only scalar changes as in (54) and (55)**. In particular, 
the maximum pressure is then proportional to 


2—20 
liz Gs 3) 
non (ae. o 
Using (15), we may express P, in physically significant quantities as 
1 
BAVA? C2 amar 
P,=| W?2 A2 x, idle ? (57) 


which is the Clemmow ballistic parameter. Note that (56) may be written in 
the form 


(58) 


By means of (58) one may interpret parameter effects with reference to either 
ballistic parameter. For example, the equivalent #=1 ballistic parameter value 
in (47) becomes *** 


(524) 
(2), a a ne) Pls (59) 


and the pressure ratio is computed according to 


_2 a 
sealer unl ales ee (69) 


It is interesting to observe that with the present formulation [7.e., equations (53) | 
of the ballistic problem, we are now able to obtain the #=0 solution from a 
standard #—=1 solution. In fact according to (59) the former is given by the 


P Ga) ale 
* Note that “‘burnt’? now becomes f3= 24) ek aeg 
5 
*xk This is a well known result due to CLEmMMow [2]. Indeed, it is CLEmmMow’s 
result which suggests the particular choices of the constants a and 6. 


xxx Here we invert the parameter ratio for convenience in the subsequent dis- 
cussions. 
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9 =1 solution with the equivalent ballistic parameter (2/P,).=0. According to 
(60) the corresponding pressure ratio is given by 


P rf sae [2 ne Ce 


The continuous dependence of the solution on the geometric parameter # in 
the range OS #1 (including zero) is thus illustrated. It is not difficult to show 
that the #=0 solution may also be obtained from the standard # = — } solution 
with the equivalent ballistic parameter value (j/F,),=0. Thus every solution 
with P,/P,=0, regardless of the value of #(—1<¥#) is, except for notation, the 
& =0 solution. 

9. Ideally progressive geometry, ® = — 1 
We are concerned here with the exceptional case #=-—1. Equations (16) 


and (17) may be written 
d PNT 6 IP. 2B ja 
eas ae 


P P Bata ? 
du 1 he (62) 
dv P\i-af P ji-a 4 
du B(3) | 2,20 | i: 
with 
y—1/ \? 
Maia Sas sal (63) 
Fee E 
The appropriate initial conditions become 
=1, 
y=0, 
biteess (5 (64) 
age 
Lf —- 
Pah 
for 4 =0. We consider a general similarity transformation 
TNE 
Ae (3) fee 
1 
(65) 


P,\™ 
=z) 
1 


where m and m are constants. The exponent 2—f in the first of equations (65) 
is chosen so as to obtain a non-critical version of equations (62)*. By use of (65) 
equation (63) may be written: 


? 


y—1 (ey V4 i 


IP Q2nr P 2 P. ue 
(z) Ba? g (66) 
d P.\i-« 
*_— p | + 0 initially. 


dp2—B 2 D—[} \I2 
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It is thus convenient to require that m—=nr. Then with n —32—2% — 1 (7.23 
m =), equations (62) and (63) become : a 
af, (Ra 
d P p27? 
M4 Ma (67) 
dv, =1(, Iz ee 
diy Puy! : 
and 
aia a ae 
Bar E 
while the transformation (65) becomes explicitly * 
es! 
Ma ioe & . en 
69 
oo oe (69) 
=(R)” 


A solution of (67) is independent of the ballistic parameter P/P, and the pres- 
sure is simply proportional to A. Thus the parameter P is fundamental to 
the limiting case #? =—1 much as the Clemmow parameter P, is fundamental 
to the case # =0. The analogy goes even farther. For if we apply the similarity 
transformation (69) to the general case of equations (16), we obtain 

dba, (By V4 
a tg P} p 


(70) 
ayy =1"(z)° 1 
d fy £, ae 
wherein 
P\i = 
p tt ro) (5) a Oua? aT 
>= 1 : (71) 
1 


Comparing (67), (68) and (70), (71) it is seen** that, except for notation, every 
solution with P/P,=0 is the ?=—1 solution. By (58), we note that P/R=0 
corresponds to P,/P,=co so that the #—=—1 solution corresponds to B/P,=0e 
while the &=0 solution corresponds to B/P,=0. We have thus covered the 


complete spectrum of ballistic parameter values, 0S a <-too, as well as the 


closed interval, —1<%<1, of geometric parameter values. 

The general discussions in §§ 4 and 5 regarding the trajectory and pressure 
characteristics carry over to the #=—1 case where many of the results become 
sharper. [See for example (30).] For = —1, the asymptotic (constant pressure) 
solution introduced in § 5 is obtained as an ordinary solution in the phase space! 
This may be seen as follows: From the general results in § 5, we infer that the 


* “Burnt” becomes w= fe 
i 

xx A similarity transformation is required here. 
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—2 
one 
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pressure ratio P/F, as a function of & might appear as in Fig. 4. In particular, 
1 


it rises to a maximum exceeding (2/)2—2% and subsequently decays to this value 
as €->co. However, as a function of x =%)é, the pressure ratio P/E curve will 
be compressed toward the initial point as %)>0, all other ballistic quantities, 
1.e., P,, being held fixed. The ee form of P/P, as a function of x is, therefore, 


the constant function P/P, = (2/y)2-2%. It is easily verified that 


1 
1 A ONo=o7 P 2 
es ME t 4 
‘¢ 2 M G 
Aarts 
— soos eee 2-20 Pt 72 
Uv M = | bt) (7 ) 


——— 


Aes = eles P*t 


is a constant pressure solution of equa- 
1 


Fig. 4. Pressure curve for #= —1 tions (6) with Pee (2/7) 2-24 and satis- 
fies the initial condition, x»=0. Since 
the piezometric efficiency for this solution is equal to one, we refer to )=—1 


as the ideally progressive geometry. We note here that (72) is an isothermal 
solution with 


It should be observed that any reasonable approximation to the #=—1 
solution requires not only that # be near —1, but according to equation (71), 
that the product (1+) (P/B)* be small. This becomes important for any con- 
siderations of a double limit ?--—1, x)~0 in an attempt to approximate 
the constant pressure solution. For with other quantities fixed, x)» 0 corresponds 
to Poo! 

10. The isothermal model and ballistic efficiencies 


Because they are often simpler to treat, many ballistic theories are founded 
on isothermal expansions. Experience (mostly experimental and computational) 
has shown that the temperature drop of the gases during burning usually entails 
but a few per cent, so that the equation of state (3) very nearly describes an iso- 
thermal expansion. However, as the kinetic energy imparted to the shot (7.e., 
the end product of the gun) is manifest in this very temperature drop, a strictly 
isothermal model would seem somewhat incongruous. In practice, the temperature 
drop is accounted for with an approximation to the temperature ratio T/T) by 
an empirical constant less than unity. This effectively reduces the force F of 
the propellant in proportion to the energy utilized. The fraction of potential 
energy of the propellant utilized by a gun is defined as the ballistic efficiency. 
Thus the validity of an isothermal solution is inextricably tied to consideration 
of the ballistic efficiency. 

CLEMMow [2] has shown that for a linear burning law, « =1, an isothermal 
expansion with a reduced force is justified wp to maximum pressure. Here we 
extend this result to all cases and show that for propellants with progressive 
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burning characteristics (7.e. #<0) an isothermal expansion with reduced force 
is justified throughout burning. 


From (21) we derive the inequality * 


f2(3)> 7) 


which is valid up to maximum pressure**. Hence using (17), (73) implies that 


(1 +9) wP—o 28 = y2 


12) 
—— 
ae ! le 2 
By Ie} 
or that 
2 ypre< (1-8) uP —8 PP (74) 
2 y 


holds up to maximum pressure. The ballistic efficiency is given by the ratio 


ee a eee 
& = GB) PO pee? (75) 


so that (74) implies that e, does not exceed y—1/y. Because of (4), the latter 
merely expresses the fact that 


(76) 


The value of y is nearly constant for all propellants varying but a few per cent 
from the nominal value 1.25. Hence isothermal solutions utilizing, say, the con- 


fe 
i 
2, 


as an approximation to ae will be universally accurate to 
0 
about 10%, up to maximum pressure. Since the constant pressure solution ob- 


tained in § 9 is an isothermal solution with 7/7)=1/y, (76) is the best possible 
result of universal character. 


stant value 


The bound (76) can be justified in general only up to maximum pressure. 
However from our results in § 5 (depicting the asymptotic pressure values for 
progressive burning propellants), we may obtain some useful extensions of (76). 
For example, the temperature ratio may be expressed as 


oe 
mp £, 
ie Ear ae me (77) 
ee ees 
One extension is obtained by generalizing (73) to 
JAN oF 
SS [aCe 
ees (z,] 2° ua 


2 
* Of course, the stronger result = (=) os +4 apples, but (73) is sufficient here. 


xx Here again we refer to the first local maximum of P. 


28* 
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which holds throughout. Then we have from (77) 


it = 1 > 


eo 
eee AL A) P 


al 

a a (79) 

throughout. Thus an isothermal solution is justified beyond maximum pressure 

in cases when the pressure P does not fall off too rapidly. Actually (21) indi- 

cates that for those cases when P approaches an asymptotic value as >, 
aa 

1.é., 0<0, the ratio* — apes = so that (77) implies that a ; as w—>oo, In 


§ Ty 


such cases we infer that (76) probably holds throughout! 
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The Equilibrium of Rods 
A. E. GREEN 


1. Introduction 


In a recent paper J. L. Er1cKSEN & C. TRUESDELL (1959) discuss the exact 
theory of stress and strain in rods and shells. They point out that stress resul- 
tants and couple resultants in shells may be regarded as derived quantities, 
defined in terms of the three-dimensional stress tensor, and that the equations 
of equilibrium for these resultants have been obtained from the three-dimensional 
equations of equilibrium for stresses by some writers. ERICKSEN & TRUESDELL 
also mention that an exact derivation of equations of equilibrium for rods based 
on the three-dimensional equations of equilibrium has not been worked out 
and the main purpose of this paper is to fill this gap. Using the explicit formulae 
for stress resultants and couples derived here, together with stress strain relations 
for large deformations, it may later be possible to obtain satisfactory expressions 
for these resultants in order to complete the theory of rods. 


2. Notation and Formulae 
We use the notation of GREEN & ZERNA (1954) and summarize the relevant 
formulae here. Consider an elastic body* which changes from an undeformed 
to a deformed position and is then maintained in equilibrium by external forces. 
Points of the deformed body are defined by a general curvilinear system of 
coordinates # where, for convenience, we use the notation 


P=f, P=n, P=e. (2.1) 


The covariant and contravariant base vectors at each point of the deformed 
body are denoted respectively by G;, G' with corresponding metric tensors 
G;;, G'’, latin indices having values 1, 2,3. If a comma denotes partial diffe- 
cn then 


G, , = I, G,, G', =— Tj, G, (22) 
where Jj’, are Christoffel symbols of the second kind and 
Mike Th, 
Liman ae Gait Sinem Gia) co 
Diyala Gre, 
Mae Ae = |Gijl- 


* The work presented here is, in fact, valid for a body which is not necessarily 
elastic. 
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The contravariant stress tensor is t’’ and the stress vector t across a surface 
whose unit normal* is wu is given by 


TT: =a, 
ps Uj sete. Te G 7 G,, (2.4) 
Vc igi 
where ,; 
UU = U; G, G;; a G; U G;. (2.5) 


If body forces are neglected ** one form of the equations of equilibrium is 
T,;,=0, Gx T=0, (2.6) 


where a comma denotes partial differentiation. From the second of equations 
(2.6) we get the symmetry property of the stress tensor. 


3. Stress resultants and couples for a rod 
We now suppose that the position vector R of any point of the deformed 
body is an analytic*** function of the coordinates €, 7 in some range of values 
of €, 7 including £=7=0, so that 


R=R*O)+EAO +4, +--, (3.1) 


where R*, A,, A, are the values of R, G,, G, respectively when £=7=0. The 
equation =7=0 defines a curve C with position vector R*(¢) which we call 
the deformed rod. The boundary S of the deformed body is defined by an equation 
of the form 


1(é, 4) =0 (3.2) 


such that €= constant are curved sections **** of this body bounded by closed 
curves. From (3.1) and (2.2) we have 


G,=A, +E A, +n IA, +>, G3) 


where ss denotes the value of Jj; when =7=0, and 4;=R*;. Without loss 
of generality we may choose the coordinates &, 7, € so that A,, A,, A; are unit 
vectors. Hence 

A, A, =A, A = Ao (3.4) 


Ay, A‘!, A denote the values of G;;, G', G respectively when £=7=0 it 
follows that 


A;,=A;-A,, Aan ae Ase As (3x5) 


and recalling (2.3) we have 


Iis1 T5330 333 = 0, (3.6) 


* The notation here for the unit normal vector differs from that in GREEN & 
ZERNA’S bcook. 


** Body forces may be included if required but for simplicity they are omitted 
here. 
fast We may assume a polynomial in € and 7 for R but this allows less freedom 
in the choice of the position vector in the undeformed body. 


Sheas If (3.1) is terminated after three terms then ¢ =constant are plane sections 
of the deformed body. 
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where Tne denotes the value of J};, when &=7=0. Also 


yr 
/G = [G, @, G,], 4 =[4, 4, Ag], (3e7) 
and we put ie r it 
H=VG/A4) =14+28y+7yI+-. (3.8) 


The force across an element of area in a ¢-surface is I, d#1 d#? so that the stress 
resultant vector m across such a surface is 
n= ff Td& dy (3-9) 


the integral being over a surface = constant bounded by the surface (3.2). 
The couple resultant vector m across the surface € = constant, about the point 
&—7—0 in the surface, is 


m = ff (R— R*) x Fdédn. (3.10) 
If dS is an element of surface area in the surface (3.2) whose unit normal 
is u then 

udS=u, G'dS = (G,d0! + G,d8) x Gd P= (@ d?— G2 dh) Gad. (3.11) 
When surface forces are prescribed on (3.2) the force on dS is tdS and using 

(2.4) and (3.11) this becomes 
tds =(de— That) ad. (3.12) 
The surface forces over the surface (3.2) may be replaced by a force p and couple q, 


measured per unit length of the rod C. An element of length of C is ds=d¢ 
so that, in view of (3.12), 


p=$ (Ldn — Tdé), 


(3.13) 
q = $(R— R*) x(Ldy — Td6), 
the integrals being taken around the curve 
6=—constanty 6f/()7)) = 0% (3.14) 
From (2.4), (3.3), (3.7) and (3.9) we obtain 

n=n A, (3.15) 

where : i ¥ 
w=VAlf Hee (i +élh + yl +---) dédn. (3.16) 


In (3.16) the double integrals are over the surface € = constant inside the surface 
(3.2). Also, from (2.4), (3-1), (3.7) and (3.10) we have 


m =A lf ox (0; en nl +...) (§A,+7A,+---)xXA,dEdn. (3.17) 
But 


A, x A; = ¢;;,.A* (3.18) 


where ¢;;,=+ V/A according as 7, 7, k is an even or odd permutation of 1, 2, 3 
and ¢,;,—=0 otherwise, so that (3.17) can be expressed in the form 


m=m,A’ (3.19) 
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where 
m, =A ff He? (6:4 ELA tliat ++) (6 tier +0 baer +) dE dN. (3-20) 
Finally, we express p and q in terms of vectors A;,, A’ respectively in the forms 


p=pA;, g=9G,4". (3.21) 


4. Equations of equilibrium for rods 
From (2.6) we have 
Sf %,,d&dn=0, Jf(R—R*)xT, ,dédn=0, (4.1) 


the integrals being over a section ¢ = constant of the deformed body bounded 
by the surface (3.2). The first of equation (4.1), together with (3.9) and (3.13), 
gives 
on 
as == (Vj Wye 
ae ty BeS0 (4.2) 
Again, from (4.1), we see that 
ra] é i pl 
Zeal (R — R*) xT dédn — i (Gs; — As) x Bedédn+ 
+ ff [((R-R*) x Bh, + (R—R*) xB} ] dé dq— 
— [f (Gx + @,xf) d&dn =0 
if we remember (3.1). Hence, using (3.9) and (3.10) this reduces to 
Se + Asxnt O (R — RY) x (dy — Tak) — ff @,xTdedn =0 


and in view of (3.13) and (2.6) we finally obtain 


om 
Os 


+A,xn+q=0. (4.3) 


Equations (4.2) and (4.3) may be written in terms of components. From 
(3.15) we have 


ae? Bere (4.4) 
where 


is n', A, Se Is A; , (4. 5) 
if we use (2.2) with =7=0. Similarly 


Om OOM; Ai 
“ig tg et (4-6) 


where t. 
13 mM, . (4. 7) 


It follows from (3.21) and (4.4)—(4.7) that equations (4.2) and (4.3) have the 
alternative forms 
ont 


1 6 i oe 
se t=O, Fen Weites He Gee 0 (4.8) 


where &2= A, ¢;=—/A and «¢,,=0 otherwise, « having the values 1, 2. 
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In the special case when the triad of unit vectors A; is orthogonal 


Ay, =A,,=Az,=0, A=1, A'=A; AA, (4.9) 


4? 


and using (2.3) and (3.6) we may write 


ae TT hae 73 

A aa ah — 133; 

oy pee eae 72 

lig a I a es EY) (4.10) 

os ph eee reany BD . 

UU ae Lo i5 

pl p12 713 

3 = 1,3 = 133 =0 


Equations (4.8) then reduce to a well-known classical form. 


5. Remarks on further developments 


In order to obtain explicit formulae for n' and m; from (3.16) and (3.20) 
we must use stress-strain relations. The position vector of any point of the 
undeformed body which corresponds to (3.1) is* 


r= 9*(C) + Ea, (0) +94, (0) +-- (5.1) 


if r is, for example, an analytic function of &, 7 in some range of values of these 
variables including =7=0. The vectors a,, a, in (5.1) depend only on € and 
are not necessarily unit vectors. The vector 


or 
as = 7 = 


(5.2) 


is not necessarily a unit vector but is tangent to the curve §=7=0 which may 
be regarded as the undeformed rod c. From (5.1) we obtain expressions for 
base vectors g;, g' in the undeformed body and hence metric tensors g;;, g’’. 
If the body is isotropic initially we then find the stress tensor t'? from formulae 
in GREEN & ZERNA. Difficulties remain in finding explicit value; for n‘, m,; and 
this question is not considered further here. However, the analysis will be similar 
to that used by Hay (1942). 
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1. Introduction 


The purpose of this article is to discuss a particular method of finding the 
solution of certain diffraction problems. Only scalar problems will be considered. 
The particular advantage of this method is that it leads quite readily to the 
exact solution as well as the asymptotic solution, with respect to large values 
of the wave number f, and that both the reflected and the diffracted wave are 
obtained simultaneously. 

The particular configurations dealt with are the circular cylinder [/], the 
sphere [2], the parabolic cylinder [3], and as a special case of the latter the 
half plane [4], [5]. Each of these problems has been dealt with in the past, 
but it is hoped that this method will throw some new light on the nature of 
these solutions. The chief advantage is that all of these problems lend them- 
selves readily to attack by this method. A complete bibliography is available 
in reference [2]. 

The problems involving convex surfaces have proved to be particularly 
obstinate and the method which was used to derive the asymptotic properties 
of the diffracted wave in all known cases was first developed by G. N. WATSON 
[6], [7]. This Watson transformation will prove to be fundamental in the present 
method. 

Another method, which has been very successful in determining the leading 
terms in the asymptotic expansions, is due to J. B. KELier [8], [9], [10]. It 
is an extension of the standard ray tracing technique, but diffracted rays are 
introduced. Unfortunately, no general proof of the validity of the method has 
as yet been furnished. The present treatment of the diffraction problems de- 
scribed above appears to have the following advantage: In the transition from 
wave optics to geometrical optics, the point (or points) on the diffracting body 
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which are relevant in the theory of ray optics automatically appear in the asymp- 
totic theory of the solutions of the rigorous diffraction problem. 

It seems that the asymptotic formulas in the Appendices are not readily 
available in the literature. 


2. Discussion of Method 


The class of problems to be discussed can be formulated as follows: « must 
satisfy the conditions 


14. 4u+Ru=o0, 
2. u+u;=0 on surface 6, 
3. wu satisfies the Sommerfeld radiation condition at oo. 


Here wu; is some incoming wave. In this article u; will always be taken to be 
a plane wave, but as will be apparent the case of a spherical wave could be 
treated in an analogous manner. Similarly the boundary condition could be 
replaced by the condition 

0 


a (™+u;)=0 on surface G. 


The solution to this problem will now be assumed to be of the form 
u=f{GH(s)ds. 
S 


Here G is the free space Green’s function, 7.e., H$)(kR) in two dimensional space 
R 


eikR 
BOE 
chosen kernel and KR is the distance from the point of observation to the surface. 
Evidently this integral satisfies the wave equation, as well as the radiation con- 
dition. /(s) must now be so selected that the boundary condition is also satisfied. 
This leads to an integral equation, but fortunately, as will be seen subsequently, 
when suitable expansion theorems exist, the solution can be found quite readily. 

/(s) is then determined in the form of a line integral or sum, which can be 
rewritten as a line or contour integral. Thus w is ultimately determined in the 
form of a multiple integral. To this multiple integral a combination of the 
method of stationary phase and saddle point integration is then applied. It 
then develops that only two points in the space over which the integration is 
performed furnish sizable contributions. One of these corresponds to the reflected 
and the other to the diffracted wave. Furthermore, on the caustic, which in 
these problems corresponds to the boundary of the shadow region, the two 
points coalesce into a single point. 

The chief difficulties arise in the evaluation of the multiple integrals. In 
order to find the points of stationary phase one must first represent the integrands 
by suitable asymptotic series. For the case of the circular cylinder and the 
sphere one has to use the Debye approximation to the Hankel functions [//]. 
For the parabolic cylinder their equivalents are derived in the appendices for 
the Weber or so-called parabolic cylinder functions. Once this has been done 
one can solve for the points of stationary phase. This will lead to several equations 
in as many unknowns, whose solution could be very cumbersome were it not 
for our knowledge of the basic laws of reflection. After this step the evaluation 


of the integrals is simple. 


in three dimensional space, © is the boundary surface, /(s) a suitably 
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3. Diffraction by a Circular Cylinder 
We shall assume that the incoming wave is a plane wave, 


“WS eikrcosé __ = ay A kr)e sages 


and the scattered wave will be represented by 


u, =f H® (kR) {(9) a9 (1) 


where 


R= |r + a — 2racos(p —9). 


Here (7, y) denote the coordinates of the point of observation, @ is the radius 
of the circular cylinder and # describes the surface of the cylinder. (See Fig. 1.) 


/ncoming plane wave I 
———— 


Fig. 1. Circular cylinder with incoming plane wave 


Use of the addition theorem [J] for the Hankel functions allows one to 
show that 


7(8) = { q” eind 


2% 4 HY (ka) * 


The resultant u, reduces to — u; for r=a. 


One can now rewrite /(#) as a contour integral: 


~ pn 0 D e—ival2 
Oy se 4 COSWU ne cos vy &e- 
@) 2m Di HI) (ka) 200 HWY (ka)siny x se (2) 
€ 
We (0) 
i 
O eine 


The contour © is shown in Fig. 2 


Thus wu, is represented as a oan integral, and one must use double stationary 
points to ene an asymptotic evaluation of the integral. To do so one must 
express H\") (ka) by means of the Debye approximation [JZ] 


HY) (k a) ww / 2 = ei n/4 et ka [cos a+ (a—m/2) sin «] 


azkacosa 
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where sin «=y/ka. Similarly one must consider the asymptotic value of H! (RR), 
which is [11] 
(1) ww |/__ 2 ilkR=a/4) 
H‘ (k R) \/ ae 


Furthermore it will be assumed throughout this section that 7 is large, so that R 
can also be approximated by 


Rwr—acos(y —#). 
Under these assumptions one can now rewrite w, in the form 


wg Sem ‘i ‘i e'8(cos a)#dd dy (3) 


—2a € 
where 


g=—hka|cos(p — 8) + cosa + (a — =) sin % +240) sina — sina. 


»-plone 


Fig. 2. Path of integration for integral (2) 


In forming g, cosy # was replaced by ~ e*”® and sinvz by — oa e-'**, These 
are the predominant terms, since as the subsequent analysis will show the region 


of interest in the v-integral is the region : ©” 0. To find points of stationary 
phase we now set nite 


eee ka[sin (vy — 8) +sina] =0, 
oe = — kal«+%—a]cosx=0. 


It is evident that the second of these two equations has two roots corresponding 
to the roots of the two factors. One now finds that there are two sets of roots 


u— Dp eR 
Phe 2a 


ut = ue 
=> +n, BO Lae 


The first of these corresponds to the reflected wave and the second to the dif- 
fracted wave. It can be observed that the former case yields the angle, which 
would be furnished by geometrical optics. In the latter case one finds that # 
is independent of the »-coordinate of the point of observation because » was 
taken to be large. One can now express g in terms of Taylor series in the 
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neighborhood of these points. Thus, to obtain the reflected wave one writes 


g=—2hasin® 4 f sin £](0- —249\)_ (4-2-8) 2(9 2EP) (q— 25%), 


Inserting this g in (3) and then evaluating the resultant Fresnel integrals asymp- 
totically for large k leads to 


an 
Sie 7. Phe) 
k(v—2 asin — 
ei ae ys (" asin 5) (4) 


In order to obtain the diffracted wave one must proceed in a somewhat more 
careful manner. The point «=2/2 is equivalent to letting y=a, but the function 


Fig. 3. Deformation of path of integration for integral (2) 


H\) (ka) which appears in the denominator of (2) has roots, as a function of » 
near that point. It can be shown that these roots are given by the expression [4 | 


; a (4m-+ 3) te 
Vea kalsa - re We =O, AO bg 


[eH a)) (yg (e ba ee 
Y—Vn 4 (R a)? 


| Ov 
One can now deform the contour of the integral (2) from G, into ©,. (See Fig. 3.) 
To do so one must justify that the contributions along 4B,CD, DE and EF 
vanish. That this is so along AB, CD and EF follows from an examination of 
the asymptotic properties of the integrand when » is large. One can show that 
the contribution along DE vanishes because on that path the integrand is an 


odd function. One can now represent f(#) as a residue summation over the 
roots of the Hankel functions: 


(ka)é eb 7/6 — im 7/2 
} (9) = a)* |Z Me ” cos Vm am é (6) 


yee * sin Divas 


and that 


Sar 7 (4m + 3) [emis (5) 
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To obtain the diffracted wave we return to (1) and insert (6). We thus obtain, 
after replacing cosy,,% by its exponential equivalent, a sum of integrals of the 
form 

ff eile cos (79) om 9] 9, 


These can be evaluated by stationary phase techniques, because the chief con- 
tribution will be due to the neighborhood of those points where 


d 

are [—cos(p — #) +8] =0. 
Furthermore, in consideration of the convergence of the series, one must take 
the range of integration from os to = for the term with + 7% and from ao 


to a for the term with —#. The resultant integral can be expressed in terms 


of Airy functions. One can then write 


Pee one 2 Beate eer ae ag) (a9) 
> Vey ie et 27m) rae (4m+3)]" rl J 


CO 

Pa 4 ,—bin/12 iky V2 2 e 

= ae e ? | tnt Ail = (4m- 3)" 
aVky Sin Vy, 2 52 (am-+3)) 


where 
[oe} 


Ai (a) = f cos(ax+ at) dx. 


0 


4. Diffraction by a Sphere 
We now turn to the problem of scattering by a sphere. The incoming wave 
will be expressed in spherical coordinates 
u, = EFF —S" (2m +1) 0" | Fa le (cos 9), 
0 
and the scattered wave can be represented by 


2m 0 


ikR 
lel £7 1(8, 9) sin ddddg (8) 


where 


IR = Ve +a*—2r acosy, 
cos y = cos # cos Hy + sin J sin Vy CoS (P — Po) - 


Here (7, %, Go) represent the spherical coordinates of the point of observation 
and (a, #, y) the spherical coordinates of the scattering sphere. 


Since the axis of symmetry of the plane wave was chosen to be the Z axis, 

gikR 
tRR 
into a series of characteristic functions of spherical waves [//]. Then, using 


one can expect that /(#, y) will be independent of gm. One can now expand 
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the orthogonality of the spherical harmonics, one can show that 


1 ak (2n+ 1) 7” F, (cos 8) 


Ba pee Ors “AY (ka) 


if 
“= —Uu Or r=a4. 


One can express /(#) as a contour integral 


Plas —ivn/2 

han 4 yP_,(cos#) e 

== : Ws a a i he 7, ) 

i) / 2n 2ni | H@(ka)cosyvx (9) 
C 


where © encloses the roots of cos yz on the positive real axis. As with the circular 
cylinder we replace H® (ka) by the Debye approximation and the Legendre 
function by [17] 


P,_;(cos 0) ~ | cap 00800 — 2/4 P=ze>o, 


~ Jy (Vr? 29) e>9Z0. 


Then (9) is put into (8), and the resultant triple integral can be evaluated by 
stationary phase techniques. The term R can be approximated by 


Rwer—acosy=r — a(cos# cos J + sin # sin Jy cos (~ — Po), 


and now the integrand contains the term e’*, whose stationary points must be 
found. Then 


g =—ka[cos# cos % + sin # sin By cos (p — Yo) + cosa + asina —(x—#) sina], 


and we have to solve the equations 


MU eer $n =9 S ie 


One root will be near «2/2, where the integrand in (9) will have poles; 
this is the term corresponding to the diffracted wave. Corresponding to the 
reflected wave are found roots at 


A ee ero ae — By 


P=; 2 2 , 
—?b a 
P=OMt2, t= a 


The first of these points corresponds to the point of geometric reflection. The 


second one can be disregarded because a= 2 %0 does not correspond to a point 


in the right half »-plane. One can now expand g in a Taylor series about the 
first point, and one obtains 


= a) 1 Onn. 
g=ha|—2asin 2 + 4 cos. sin By(p — go)? + 


— sin %o [9 ee 3 [x = + 3) 


No Oe Tt By \2 sls MARTY 
+- gees 3 (9% _— <2) => sin (= 


2 2 2 
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The evaluation of (8) now leads to the expression for the asymptotic form of 

the reflected wave, namely 

a ih(r—2a sins) 

Us R 2R Oo’, (10) 
To obtain the diffracted wave one can return to (9) and integrate over the 

residues due to the roots of the Hankel functions in the denominator. This is 

done in a manner completely analogous to the process carried out in the problem 

involving the circular cylinder. The result (9) can then be rewritten in the form 


x DOE a Pe nee ek — g/4) etm m4 
FG ea) é » (—)" 0m? COS(Ym— 7/4) e (11) 


V2zsin 0 - fe Coe esy oe 


Here the y,, denote the roots of the Hankel functions. (11) is then inserted in 
(8) and in order to evaluate the resultant integrals one can split the cosine term 
into its exponential components. Only the term e '’**'7/4 contributes points 
of stationary phase, so that the other term may be neglected. In order to evaluate 
(11) one must evaluate integrals of the form 


20 7% 


f fe tbacesy+om 9) Vsin Fad dg. (12) 
0 0 
The points of stationary phase are furnished by the roots of the equations 
$5 [cosy + 3] =0 


¢G anh —— 
bee [cosy + 8] =0 


Pe Lor j=0,— == 
which are ieee The phase must be chosen properly, to allow 
P=Potm, 0=0,— 2 


the series to converge. The q-integral then becomes a Fresnel integral and the 
#-integral an Airy integral. Actually the g-integration can be carried out exactly 
in terms of Bessel functions, which are then approximated by their asymptotic 
series. It appears that little is gained in performing the exact integration. The 
above process is valid only if 0<eS%)Sa—e. For Jj)=0 or a, there is no 
point of stationary phase in gy. Upon performing the resultant integrations one 
obtains the diffracted wave. 


co 


RY A —15n/12 —47 tym (n—O + 1m (m+9, 
—e'*"(2ka)se 15 21/ » (—)" vm? [e Vn ( )—j etm | +90) 


krVsin 8 COS Vy 20 ee (4m+ 3) 


0 


x eed Se (4m + 3) | ; 
Fe per anc) 
ane (16h abe *7"8 » (=)" vibe! 008 (98/4) 
kh y Vsin De COS Vy, TU 2 (4 m-+3)|" 


se 2 — } f al 
S< a 8 (4m +3)| { 6 
Oa Oat kc 
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To find the diffracted wave for J) =0 we return to (12), where the @ integration 
is performed directly since now the integrand is independent of m. The #-integral 
is evaluated in the same manner as before, and we obtain 


[o.e) 


& akyr B  —15m/12 __\M, —k tm ° 
Tae Mi (2ka)ie Duin: a a Ail pe yes 3)". (14) 
a 608% 0| 2 (4m-+3)|" 


The case 0)=2 yields a similar expression. To evaluate it one must return 
to (41), but the complex term e‘’°~7/* must be used now, because the other 
term provides no point of stationary phase. 


5. Diffraction by a Parabolic Cylinder 
The technique applied in the earlier sections can be applied to this problem 
as well. But now one must operate with the appropriate functions and expansions 
in parabolic coordinates. It can be shown [/2] that 


—a+1co 


V2e-%*7 14 =< ‘ 1 1 1 ] 
(1) = a wen Are Bee eis 
” eh) I ant ee | r| areas lea. 
x EO. (2V =the) D3. 1 2) x (15) 
x EP_,(2V—tky.) Dy, 4(2V—ikn,) ds 
and 
1 jh pcome tee 1 1 
tk(xcosat+ysina] _ ee Pe! 1 1 ae) ew 
i V2 sin « 2m 4 o- | r| oa 4 | 2) zs | 4+ 7) x 
¥ Pes AP oe ae (16) 
x (—}{tan&)\ "£9, (2 V—7kL) EO_,(2V—ikn) ds 
2 2 a / 
where 


OSasm, \lz0, oS|yil<e, |o|<f+ 5. 


The parabolic coordinates are defined by 


x=C—y, y= H2yln 
and in (15) the symbols ¢_, ¢. are to be interpreted as the ¢ coordinates of two 
points where ¢_ is the smaller and ¢. the larger, and similarly for 7 Lo) Sage teas 
the distance between the two points. The function D,(Z) is a solution of the 
differential equation 


yt 


2 
oe 3 —Fly=0 


and can be defined in terms of the Whittaker function 
ye (Z?\-3 z 
D,(Z) = 2" (=) Mort 21(5)- 


The functions E£{)(Z) are solutions of the same equation and are defined by 
1+-» 

(2) Pe 1—l—» 

E? (Z) 2 F(t = ){D,(Z) + (—)'D,(—Z)},  1=0,1. 
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It can be seen from this definition that E!(Z) is an even function of Z, and 
E (Z) is odd. The E and D functions play the same role in parabolic coordinates 
that the Bessel and Hankel functions play in polar coordinates, respectively. 
In other words, one associates standing waves with the FE functions and cylindrical 
waves with the D functions. This may be observed in (15) and (16). 

The scattering parabolic cylinder 
will be defined by €=€, so that the 
scattered wave will be representable 

s (see Fig. 4) 


EE, Te 


“= J HY (k R) f(y) d V—27k 7. (17) 


Here & is the distance from the 
point of observation to the cylinder. 

Again f(y) must be so selected 
that on the surface uw, cancels out 
the incoming plane wave. This leads 
again to an integral equation for / (7), 
which can be solved without too much difficulty by assuming a suitable form. 
Thus it appears reasonable to suppose that 


Direction of plane wave 


Fig. 4. Parabolic cylinder with incoming plane wave 


—o+i00 


A,(s) (ton 3) Fah ,(2V—ikn) 
/(n) se ; 


Ds 2s—3 (2V—a Ra) 


PS (18) 


—a—1ico 


The E functions in the numerator are orthogonal, since one is even and the other 
odd, and in the denominator we require a D function, which will cancel the 
corresponding D function in the numerator of (15), so that when evaluating 
(17) for a point of observation on the surface the resultant expression will resemble 
the plane wave (16). 

To determine the unknown terms 4;(s) in (48) one inserts (18) and (15) 
n (47), and one has to use the integral formulas 


Wee hae hetero iy 
0 
Van 2+ ED), (2V—th 1) 
—j = 0,1. 
4(s 0) P| s4 a 


Here mp is the 7 coordinate of the point of observation and 


0, %o) 


No» °°). 


Ne=N, Ns =No | 


for the interval of integration ( 
em te ales aren 


These formulas can easily be derived from general integration formulas given 
in the work by Bucuuo1z [12]. 


29* 
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(17) then reduces to a double line integral 


—a+ioo —a+100 


Biers sien | ‘i r(s+4 fii Apo Or ,(2V—tk£) x 


paasi 2) ayer —W2a2h) po (9 iky)dsde. 
m Es as AG) [tan | 27 1 ( V No) 


One can now close the contour of integration for the variable s, using the appro- 
priate asymptotic formulas [72], thereby evaluating the integral by means of 
its residue at the one singualr point s=7t. One is then left, for €>=¢,, with 


—a+100 


wa BEVEL f T+ eta) Gl=ie) x 


—a—100 


x Ea, (2V= ihm) Ai(s) [tan Z) aba. 


Since this is the value of the scattered wave on the surface, it must reduce to 
(16), and one obtains by comparison 


l ae 
A ir stot 5)ai 
A,(s) = a “Vsin a 25 
so that finally 
mete | 1 | l a \~?s (2) 25 bw 
eee gel r nae 2) (tn 3] Evsas—a\ Vik n) 
f(y) Pipe : i as. = (49) 
“Vsin x Es = o Dag 2 Tec.) 


—o-10 


Furthermore, one can show, using the linear dependence between the & and 
the D functions, that 


—a+tioo 


(es 


"4 4 a\—2s Dor (2V=2 k ”) 
ee lhe ae : BSL ay Fao: 
ii) = 5 i r(4 —2s)(tan ©) amen Sree 


—o—100 


The next step is to combine (17) and (20) and to seek the stationary points 
in that double integral. To do so one must express the D functions in terms of 
appropriate asymptotic formulas, which are developed in the appendices. The 
details of this analysis are fundamentally the same as the corresponding steps 
in the analysis of the sphere and circular cylinder, but become much clumsier 
in the case of the reflected wave. Again two pairs of roots are found for the 
equations yielding the points of stationary phase. One pair corresponds to the 
reflected wave and the other to the diffracted wave. As was to be expected, 
the points yielding the diffracted wave are given by 


oR : 
Bi tat ee eae 


Here the parameter ¢ represents the arclength of the parabola, and the important 
point is the point where the tangent ray from the point of observation to the 
parabola touches the surface. In the s integral we are again close to the residues 
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of the integral. The reflected wave depends on the points where 


Ohaaas “—t Ee: 
Sa eae 28=1kC,|1= 


8 being the angular coordinate of the point of observation. The border of the 
shadow region and the illuminated region corresponds to .=#, where the reflected 
and diffracted wave merge into one. By a direct but lengthy calculation one 
now finds that the reflected wave is given by 


at & 


a 0 ee 2 sin (a—#®) sin } 
Konan 4c y si a+ d i, até o 
pros CSC - DSece sin } 
2 2 
x exptk ie 0 sin 8 [cse oe — sin *4*) +osin?«sin? #) (21) 


where @ is the distance of the point of reflection from the origin. 

The diffracted wave can be obtained much more easily. To do so we return 
to (20) and find by a direct expansion over the residues of the D function in 
the denominator using the asymptotic formulas (A 8), (A 9), (A 13) that 


m & Son 2 Sm 
Vane (et,)8 yO) (ee ale 
i) 16 sin a ie $47. , thnk es 
: ce Be ) sn( 1+ 2s ‘| 
m 
7) (Uesarerl ly ace) ee 
ae 4 a (2+ fe z)| (22) 


This expression is now inserted in (17), and the Hankel function is replaced 
by its asymptotic form. If that integral is evaluated by the method of stationary 
phase one finds that one is concerned with the point where the function 


Ga Bart 
is stationary. It follows that 


og ai . at 
Sa le 


dyn \ at eyn 
vanishes at the point where oe = — 1, that is, the point of tangency. Furthermore 
Ohgh = WO ue ea 
CV = GWE (nes)! 


at that point. Again one is concerned with Airy integrals, and one obtains 
finally 


——IUe 


int 0 
0 
Gs, Die 256 sina (k £q)8 we ae ae 
s “itl Sata |/ EIB = 


1 


(4m- eh (14 bn) (1+ 4 


(4m -+3)]'t. 


— 
p— 
iW) 
iS 
ses 


8 


434 Harry HocuHstaptT: 


The subscripts 0 refer to the point where the tangent ray from the point of 
observation to the parabola touches the surface, and g, is the distance of that 


point from the origin. 
6. Diffraction by a Half Plane 


The solution to the problem of diffraction by a half plane can be obtained 
by returning to equations (17) and (20) and letting ¢,+0. Then, by use of 


— feats 
D_35-4(0) = T(is+3)’ 
(17) is transformed into 
—o+io . 
1 3 s+4 a@\—2s 4 Noa 
— =e ait = IN = Decne tk as 
i(n) 2nVasin « [ i S| (tan 3) (5 s 2 3 ( | n) 


This integral can now be evaluated in terms of the residues of the integrand 
in the left half plane, so that 


es ach ae Oe eer Neve a 
1 (n) anee N! 2 a (tan 5) Ri sate N) D_ an—2 (2 V —ikn); 


and, using the integral representation [J/] 


(=e t 
—2)—tkynt— — 2N 

2 WL 
3 dt, 


P'(2+2N) D_sy-»(2¥—iky) =e" [e 


0 


one can sum the resultant exponential series to produce the closed form 


Hee —i(tan %) wnt == <i t— + sec : tat. 


V2asina 


This integral can be evaluated directly, by means of the previous integral repre- 
sentation, so that, finally 


isin eee ee 
Hn) = ce D_,(2 cos $Y — ikn). (24) 
Now one can rewrite (17) in the form 
u, = Baie {i HH (k R)e oan -2(2 cos > V—iky n) 4 \—2tkn. (25) 
2Va J 


This problem has been solved by similar techniques [13], [14], and in particular 
a very similar result was found in reference [13]. Well known asymptotic formulas 


can now be derived from (25). In particular the following representation [11] 
is helpful 


D_,(2)=e* {22 i Oe he de 
2/y2 


This can be represented in terms of the standard Fresnel functions. 
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Appendix I 


The function D_»,_1 ,(2V—tke, ) has, as a function of s, roots in the left 
half plane. These can be located by finding a suitable asymptotic form for large k, 
which will be done by LancEr’s method [15], [16]. The function D_»,_, (2//—ik¢,) 
satisfies the differential equation 


a 2) 
=e (2s + 7) 2 =o: (A1) 


New dependent and independent variables are now introduced by the sub- 
stitutions 


Z=2)—2sy(0), y= Ve y=ypu, 
leading to the equation 
d? 
Tex + 165% (1 — ge) (pu =0. (A2) 


g, which as yet is undefined, is now determined by 


In this manner the equation is reduced to Airy’s differential equation 


au 


de + 16s?¢u=0. 


One now selects the solution that tends to 0 as (ov, namely 
u=L Aif[—f(— 16s?)8] (A 4) 


where L is a suitable constant, independent of ¢. 
co 
TT (9) £f cos (= B+ ax\dx, 
m 3 
0 
and for large negative arguments 


It is also known that for large z 


De 2 Ven ia wails 
—_ sS=—% a 


From (A 3) one now finds that 


° aie te 
2a cher < py? —1 —cosh4 
20 [Ve—1dp=* 3 
i 


so that for p large 
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If one now assumes that » is large, and combines (A 4), (A 5) and (A 6), one 
can evaluate L, and one finds that 


5 I — he os fas roan 
Lata Te (s—}) Ins li i 


also one finds, for pa 1, that 


th 1 in Ss ns—sln2+anis 
Biro (A 7) 


One can now locate the roots of this function, and one finds 


| 2 (4 m-+3) | 
gras! Sab ae eee eel, M=O0, Ie 252%. (A 8) 
and 
ta) ae — 93 2 M12 5 (Sm 4) Insm SmIn2+ias 
my (2) ak a) ee ee : y, 
OSs @ a( V YZ oN m (kE,)8 
a yn |32 (4m 4 3))"- (A 9) 


Far away from the roots the approximation takes the form 


xem RO ee in/4 —sIns+ins 
D4 2) —th) 
. V2 2° one 


(A 10) 


Appendix II 


One can analyze the function D,,_, 4(2)/ —ikn) in a similar manner as was 
used in Appendix I. Dz,—4(2\/—ikn) satisfies the differential equation 


ay EE 
ae ee 7) 9=0. (A 11) 
New variables are now introduced by the substitutions 


Z=2V—2590, Wea ee y= Pu, 


and ¢ is defined by 


De ead 
a laa 


2 ey oVite+sinh g 
3 ® 


=) 


so that 
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Thus (A 11) is reduced to the Airy differential equation 


du 


One can again combine (A 4), (A 5) and (A 6) to show that for large k 


ler e2sS—sins+nis+ni/4 


oe. ee oat 
Pr; 25) Dac4(2) a) 2 (ik +25)! 


x exp {ikn (hn + 25) +2sin([/22 Yt 2a}. 


(A 12) 


DS 


In particular one can now show the heavy dependence of (A 12) on the arc- 
length ¢ to the point (¢,,7) on the parabola from the origin, when S lies near 
the points (A 8). One can show that the arclength is given by 


af JorSer=taerrtae sm [f+ YZ 


so that (A 12) ees to 


= 42 Sm—Sm IN Sm +11 Sm +1 1/4 


ip 
P= 254) Da 4 (2 V—ikn) ee a= a 
DENG Ais n+25,,)4 


x expihk {i+ One, in( | + +2 AV, (A 13) 


where 
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Representation Theorems for Solutions 
of Linear Partial Differential Equations 
in Three Variables 


JOSEPHINE MITCHELL 


Communicated by S. BERGMAN 


§ 1. Introduction 
BERGMAN has discussed the properties of integral operators which generate 
solutions of linear partial differential equations with analytic coefficients in the case 
of 2 and 3 real variables [/—6]; see also KrEyszic [8, 9]. In the case of 3 variables 
they consider the equation 


(4) Ayp+F(r*)y=0, 


where F is an entire function of 7?=x?-+ y?4 22, X=(x, y, z) a point in three 
dimensional 4%, y, z-space E,— x, vy, z realt. The methods and results obtained 
in the theory of harmonic functions of 3 variables can be utilized in the study 
of solutions of (1), but first the theory of harmonic functions of three variables 
had to be developed as a basis for the theory of functions satisfying (1). See 
[1, 2,3]. These results on harmonic functions are of interest in themselves. 


The Whittaker-Bergman operator ([Z], [5] p. 467): 
(2) H(X) = BUf, 8, X_) = (1/2008) f fw, 0) ate, 


transforms an analytic function / (called the B,-associate? of H) of two complex 
variables w and ¢ into a harmonic function of three variables for X in the neigh- 
borhood of a point X,. Here 

u(X,6) =*+Z2042*C4, 

Saal’) > Sa OY =) 

and & is a closed rectifiable curve in the ¢-plane. The representation (2) at a 
point X, where / is regular is valid at first in the small (that is, in a sufficiently 
small neighborhood of the point X9). BERGMAN has shown how we can use the 


correspondence /(w,¢)—>A(X) for the study of harmonic functions in the large, 
in the case where f{(w,¢) is a rational or algebraic function or the integral of 


U 
(3) ‘ 


1 Recently BERGMAN has investigated also equations Ay + F(y, z)y=0, where F 


is an entire function of y, z. 
2 We note that the B,-associate of a given harmonic function depends upon the 


choice of the coordinate system. See [5], p. 467. 
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an algebraic function. The harmonic function defined by (2) is not regular 
throughout £,. For example, if / is a rational function, the space £, is sub- 
divided into a finite number of regions (called domains of associations) which 
depend upon f and the curve & in the ¢-plane, such that within each region (2) 
defines a regular harmonic function. As the point X moves from one region 
to another, singular points of / enter or leave the interior of 2, and a new harmonic 
function is defined. . 
One problem in the study of integral operators is to obtain integral relations 
which generalize the residue and Abel theorems in the theory of analytic functions 


of one complex variable [5]. In this case we consider harmonic vectors H= 
(H,, H,, Hz), that is, vector fields in three-dimensional space which satisfy the 
system of equations # 


(4) divH=0, curlH=0 
(see [5], p. 488). In particular we consider harmonic vectors, 
5) H(X) = BU, 2, X), 
whose first component H, is given by (2) and# 
B= ae fi fuc) 4 ac+ey 
(6) ito = ty [1 ee « 
iret 
(7) Aiton prea 
st 


where 9 is a rectifiable curve lying in the regularity domain of H with initial 
point X9 and end-point X. In case f is a rational B;-associate: 


(u, C) 
160) = oe ty 


P and Q being polynomials in w and ¢, BERGMAN [5] has shown that there exist 
simple relations between a linear combination of the transcendental functions 
Ay(X, X9) and a linear combination of certain integrals of algebraic functions 


of one complex variable, 
ff . PO) dt 
{leo O)/ae ‘), ime 


representing a different class of transcendental functions. In these relations 
the arguments x, y, z of H and argument of the function of one complex variable 
are connected by the algebraic equation 


(8) Qlu(X, 0), ¢] =o. 


3 It should be noted that these vectors represent only a part of the class of vectors 
whose three components are harmonic. 

* We note that the components H, and H, to a given H, are determined only 
within Re[g(y+7z)] and —Im [g(y +22)], respectively. Here g is an analytic function 
of the complex variable y + iz, regular at the origin. 
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In this paper these relations between integrals are generalized to the case 
where / is an algebraic function. In § 2 a relation is obtained between a linear 
combination of certain transcendental functions increased by period functions 
(ef. (2.24) ]° and integrals of certain algebraic functions of one complex variable, 
also, increased by period functions (2.27). In §3 examples are given to illustrate 
the general theory. In § 4 an analogous result is obtained in the case of rational 
B;-associates for solutions of equation (1). 

Notation: 

In the following, bold-faced Italic characters denote vectors, while German 
special characters denote manifolds. Superscripts on manifolds lying in EF, 


_ denote their real dimensions. The notation E[...] means the set of points satisfy- 
ing the condition in the brackets. 


§ 2. Integral theorem for harmonic vectors with algebraic associates 

1. The domain of regularity of a function with an algebraic Bs3-associate. The 
first step in our investigations consists in determining the domain where certain 
harmonic functions are defined by (1.2), in the case of an algebraic B,-associate / 
(see Theorem 2.1). Let f be an algebraic function 

_ p(u,¢,S) 

(1) (6,8) = PES), 
where # is a polynomial in uw, € and S andg in wand¢; u,¢ and S are connected 
by an irreducible algebraic equation 


(2) v(u,C, S) = 2, p(w, 2) S*=0; 
v=0 
p, are polynomials in uw and ¢. Under the transformation (1.3), v goes into 
(3) u(u, 6, S) =f-*V(C, X, S) =C7 Digs, X) S” 
v=0 


(here #y is a positive integer and g, are polynomials in ¢ and X). For fixed X, 
V(¢, X, S)=0, considered as an equation in S defines for ¢,==0 a Riemann 
surface t(X) of m sheets and of genus @(X), lying over the ¢-plane. The genus 
o(X)=e is in general (that is, except for some values of X noted below) inde- 
pendent of X. The branch points e;(X) of }i(X) which lhe in the finite ¢-plane 
are included among the roots of g,(¢, X)=0 and roots of the discriminant 
d(c¢, X)=0 of (3). We assume d(C, X)==0. The topological structure of W(X) 
is different from that described above in the following cases: 

i. If for any X, q,,(¢, X) is zero for all €, then R(X) has less than sheets. 
This situation occurs if g, has a factor, say b(X), independent of €. Since 0(X) 
is a complex function, )(X)=0 represents in general two equations Re 6(X)=0, 
Im 0(X)=0 and hence defines a curve in £3. Since qg,(¢, X) is a polynomial 
in € and X, it can have only a finite number of factors independent of ¢, such 
as b(X). Hence there are at most a finite number of curves in £3 such that for 
X on these curves, the surface {t(X) has less than m sheets. Let Sj be the set 
of points in E; lying in these curves. 


5 (2.21) means Formula (21) of § 2. 


442 JOSEPHINE MITCHELL: 


ii. The discriminant d(¢, X) is obtained by eliminating S from HANGS ASS A ea) 
and V;(¢, S, X)=0 and hence is an entire rational function in € and X. Thus 
there are at most a finite number of curves in £; along which d(C, X) is zero 
for all. Along such curves V(¢, S, X) has a repeated irreducible factor involving 
S, cf. [7], and the structure of }t(X) changes. Call ©} this set of points. 

iii. Two branch points coincide if g,(¢,X)=0 or d(¢, X)=0 have multiple 
roots, that is, two branch points coincide along the set of points X for which 
dn(C, X)=0 and 4q,(¢, X)/eC=0 or d(C, X)=0 and d(C, X)/0C=0. Such a set 
of points (G4) usually lies on curves in £; for if¢ is eliminated from, say, ¢,,(¢, X) =0 
and 0q,(¢, X)/e¢=0, the resulting function is in general a complex function 
of X and hence defines a curve. Let 


(4) Gi = E[X|XES} a Spores) 
Under transformation (1.3) suppose that f(u,¢, S) goes into 
(GS fee) 
EGe ote) ee 
where 
(6) QO(C, X) =CPq(u,¢), 


p, an integer. The function F(Z, S, X) is determined in a Riemann domain 3, 
that is, in a domain consisting of m space sheets. If we denote by 3(¢) the 
Riemann domain of the function F(¢, S, X), ¢ fixed, then for every X except 
for X belonging to G!= Givi Sj GS}, the domain has the same topological 
structure and therefore for X €R*— G! the genus 0,(6)=o0 of RF(C) is constant. 

Let & be a rectifiable curve in the ¢-plane, not going through the origin, 
with initial and end-points ¢, and ¢, respectively and consider® 


(7) f F(t, S, X) 4. 


We determine the domain of regularity of (7). For fixed €--0 the equation 
Q(¢, X)=0 or g(u,¢)=0, which may be written in the form 

ag(0) [] [w—u,(0)]=0, w=ulé, X) 
lefs(43) 5-0 


that is 


defines a finite set of straight lines JU (C)"In72.) For ERS YA(C) trace a finite 
set of surfaces 
(8) qe aie M2 (6) = B[X] OG, X) =0, CES). 


For X €X? the integrand of (7) is undefined. 


ae note that by changing & the function given by (7) can be analytically 
continued. One of the problems in our approach is to study the properties of the 
function obtained by analytic continuation. 
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Similarly for fixed ¢ the equations g,(¢,X)=0 and d(f, X)=0 define in 


general curves Bj (¢) and %}(£) respectively and as ¢ traces 2, they trace sur- 
faces in Ey: 


(9) B= E[X|q,(¢,X)=0 or dll, X)=0, CER]. 


For C€%, O(¢, X)=0 only at those points X which by (8) belong to Q?. 
For other values of X ¢N?U G1, F(C, S, X) is a regular function of € and hence 
(7) exists at such points. 

Hence follows 


Theorem 2,1. Let F(¢, S, X) be the rational function defined by (5); €,S 
and X being connected by V(C, S, X)=0. Then the integral (7) represents a function 
which is defined on the Riemann domain 8 at all points except when X CYISA. 

2. Behavior in the large of functions H,(X) with algebraic associates. According 
to a result due to WEIERSTRASS [10], fy =F/C can be expressed for all X € 2 — SG! 


in the form 
v(X) 


BESS DGC, 5:6, 5; K+ a SR, 5, X) — 
(10) es rr 
— YU(X) Gale, S, X) — fel) GalG, 8, XY], 


where (¢,, S,) are the positions of the poles of Ay (¢,, S,, x) _ ¢ on the Riemann 


surface R(X) (¢ is the local uniformizing variable on K(X ie c,(X) is the residue 
of Fy at (¢,, S,) and G,, G,, G are Weierstrass integrands of first, second and 
third kind respectively. The functions c,(X), f,(X) and f,(X) are algebraic since 
they are obtained by algebraic operations from algebraic functions; R,(¢, S, X) 
is a rational function of €, S and X, 7.e., it is an algebraic function of € and X. 
Pet 

2Wug(X) , 20,8 (X) 
(11) 2Mug(X), 21a p(X) 

Lacy, S35 ka)s 23(6,, See 76) 


(a, B=1,...,0; v=1,...,7/(X)) represent systems of primitive periods of Abelian 
integrals of first, second and third kind respectively. If c,(X) and cg(X) are 
the cuts and conjugate cuts introduced on }#t(X) to transform it into a simply 
connected surface t’(X), then 


2. g(X) = ile,(X) eX), 2Wa p(X) ear (en( eX") 
(12) 2%. (X) =I, CARSON 273.8 (X) el (c, (X), X) 
MgC. Sp, X) =T(cg(X)3o,,5,3X),  Qall,, SX) =1(c(X),¢,, S,; X), 
where 
TC GAGS, x) ae 
6 
(13) H,(C, X) = J Ga(6, S, X) de 


TG OL) GG, 55, 0,4.) ae 
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for C an open curve joining two points of R(X) and not going through a 
branch point. rie 

Now let & be a closed regular curve on 3i(X) which does not meet the origin, 
the branch points of 9(X) or the poles of F(¢, S, X). Assuming that & inter- 
sects the cuts c,(X), cg(X) only a finite number of times and at each intersection 
crosses the cut (cf. [7]), we define m,(X) as the difference between the number 
of times & crosses cg(X) from the negative to the positive side and the number 
of times it crosses in the opposite direction; similarly define an integer 1,(X) 
with respect to c,(X). Then the Abelian integrals of the first and second kind 
over & give 


TS re) = [2 (X) Wp (X) — 2n3(X) p(X) ] = Sales) 
(14) : 
Ti (@,X) = ¥) [2mg(X) map (X) + 25 (X) mhp(X)] = See(X)- 


pal 
If loops J,(X) are introduced around the points (¢,, S,) on the surface }t’(X) 
to remove the points at which the integrand /(¢, S, X) has poles, the resulting 
surface §t/’(X) is also simply-connected (cf. [7]). If m;’(X) is the difference in 
the number of times & intersects /,(X) from the negative to the positive side 
and vice versa, then the integral of the third kind gives 


Q joe 
+2nin, (X)=S,(,, S,; X) 


(v=1,...,7(X)). Thus from (10), (13), (14) and (15) 


dt r (X) eee pes oy ’ 
[ FG. S,X) = % [eh (X) 5,6, 5X) +R] — 


— ¥ TEX) Sae(¥) = fa(8) Sio(¥1, 


where c,(X)=c,(X) if (¢,, S,) lies inside 2 and 0 otherwise and R,(X) is an 
algebraic functions of X 7. Consequently 


Theorem 2.2. An integral of the form (7) with an algebraic associate F(£, S, X) 
can be represented as a linear combination with algebraic coefficients of normal 
integrals of the first, second and third kinds associated with the Riemann surface 
N(X) of the algebraic equation V(C, S, X)=0, increased by an algebraic function, 
as given by (16). 

The expression (16) is in general a transcendental function of X defined in &3. 
(2? is a domain with infinitely many space sheets defined over the Riemann 
domain }%.) BerGMAN [2] characterized these functions by showing that all 
the period functions belonging to a given algebraic equation, considered as a function 
of x (or y or z), satisfy a system of ordinary differential equations with algebraic 
coefficients and can be singular only if two branch points coincide, that is, only 
along the set G'. It is possible to cut 28 by a piece of a surface Z? such that in 
the domain 23= ht — J? every closed curve can be reduced to a point. In this 


’ BERGMAN Called the functions (16) ‘period functions’’. 
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case the period functions w,4(X), @,,(X), 7, a(X) and »,,(X) are single-valued 
in D*. If, however, starting at X) and moving along a closed curve which inter- 
sects 2? and which cannot be reduced continuously to a point without inter- 
secting ©1, we return to Xy, each of the period functions will assume a new 
value. Similarly for Q,(¢,, S,; X). Following BERGMAN [2] we may say that the 
totality of functions represented by integral (7) with common algebraic associate 
F(C, S, X) belong to the same category. 


3. Expression for the integral of a harmonic vector with an algebraic associate. 
Let 1 be a closed sufficiently smooth oriented curve lying in 9, — S}, subject 
to eee restrictions formulated in Paragraph 4: 


(17) Sn EN IG) ys \ ee (5) Os S = Sol, 
and consider the integral 
> =e 
I,=2ni | H(X)-dX 
= fix {FCS 4) F+ pity fC+oF Hae ica C) F 
aid a 

In general $1 will meet the surfaces 9? and Y? introduced in (8) and (9) in a set 
of points X, (supposed finite) (k=O, » 8), thus subdividing the curve 9! into 


a finite number, s, of segments 9%, G!= oa %; With end-points X,_,, X, such 


that for X in the interior of $4, & does ae intersect any of the branch points 
e; (X) and the number of poles of F(¢, S, X) inside & is constant. The branch 
points of R(X) are defined by q,(¢, X)=0 or d(¢, X)=0. If a branch point 
also lies on &, then X € BY? by definition of B?. Hence X is one of the end-points 
X, of the segments 9;. Similarly if X is in the interior of 3%, no pole of F(¢, S, X) 
meets &%, that is, no zero (;(X) of Q(¢, X) meets &. Also since the zeros ¢,(X) 
vary continuously as ¢ moves continuously on % and @ is a closed curve on }i(X), 
¢;(X) can move in or out of & only for a value of X for which ¢;(X) lies on &. 
Thus for X in the interior of 9; there exists a sequence 7;,() of the zeros of 
Q(¢, X) such that ¢,(X), »€ 7, (&), lie inside 2 and the remaining zeros lie out- 
side &. 


Thus we consider 
s 


L= NOU Xe 

(19) = DA Xe Xia) 

where u 

(20) A (X,,, AG 4) 2307 si H)(X)- dX. 
Si 


By means of (16) 


(21 = 31 Bo. of [Selle So: 2200) ak +8, (x) -aX] - 
= 3) J [Ss(X) fAX) AX — S200) F120) -aX]} 
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where ¢,, e) Fa and R, are vectors whose first components are given by the corre- 
sponding functions in (16). 

4. Relation between the functions A(X, X,) and integrals of certain transcendental 
functions of a complex variable. For fixed ¢ the irreducible algebraic equation (2) 
defines an m-sheeted Riemann surface §t,(¢) of genus @, lying in the w-plane. 
The finite branch points 0,(¢) of 3, (¢) are included in the zeros of the leading 
coefficient ~,(u,C) of (2) or of the discriminant d,(u,¢) of (2). As in the case of 
the Riemann surface $t(X) the number of sheets is less than m if for some ¢ the 
polynomial #,,(w,¢) is zero for all u, which could happen for only a finite set 
S{ of points ¢. Similarly the number of branch points or their multiplicity and 
hence the genus of 3,(¢) could change for only a finite set Sj in the ¢-plane. 
Assume that the curve & intoduced above does not meet S;u S3. 


Interchanging the order of integration in (18) yields 


(22) ie: [+ [fs S) du 

g 3(C) 
where du=dx+4idy(€+0)+$dz(€ —C74) and 9(¢) is the projection of 31 on 
the u-plane: 


(23) O(C) = Elulu=u(X,¢), XES'). 


We assume that .§1 in (17) can be chosen so that §(¢) bounds a domain on ‘i (¢) 
for all €€ & except possibly a finite number. If we set €=re'", then for r+1, 


(Z) is the curve eeae! ie Im a cme (y2 422), 

lv ue 

Since the polynomial #,,(u, ¢) of (2) fon g,(X, C) of (3) divided by a certain 

power ¢! of ¢ and with uw replaced by u(X,¢), f, [w(X, ¢),€]=0 and q,(X, 6)=0 

define the same set of points Bi(¢) in E;. Also if S,(w,¢) is a root of (2), then 
Si (X,O—Si[u (X50), C] is arreot of V(C,.S).)ee0 


O=v[u,C, S,(u,C)] =v[u(X, 2), 0, S,(w(X, 2), ESN 4 IRA Y Re. 


since v(u,¢, S)=v[u(X,¢),C, S]=C-”V(C, X, S) and hence d(X,2)=0 and 
d,[u(X,¢),¢]=0 define the same set of points %3(¢) in E;. Let [€,] be the set 
of points on & which correspond to the set X, in which 9’ meets 22 and B?; 
the set [¢;] is defined as follows: X,,€ B2U N? implies X,€ BT (f)u BZ (Oho M2 (0) i 
some €=C;€& [cf. (8) and (9)]. Since Q is an oriented curve, the ¢; can be enumer- 
ated in the order in which we meet them moving continuously along 2, and & 
will be subdivided into say s’ parts having end-points €;_, and €;, fy=Cy. 
According to the Weierstrass formula used in (10) 


f(u,¢, X) = Va, (us. Sy, S342 D Re, 8,0) — 
(24) a 
= 7 [ga (2) Gy, (u, S, ¢) ae BAC) Gia(u, S, ail 
where (u,, S,) are the position of the poles of f (uy, Sy, 6) a t’ (t’ is local uniform- 


izing variable on 9, (¢)) and the functions d,, G,, G,,, Gia, &,, 2, have similar inter- 
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pretations to the corresponding functions in (10). Also R,,,d,,g, and gi are 
algebraic functions of £. By a proof similar to that for formula (21) the integral J, 
of (22) can be expressed in the form 


=>} » f a0 by, $0) +R,,(¢)] S — 
j=1 (ve F (9) 
(25) b 
J (ehl€) Toe) — galt) tO] 1, 
a=1 
where R, , (€), g,,(¢) and gi, (¢) are algebraic functions of C; TF, = AG et POA INR A Sa 


cae 
are transcendental functions corresponding to the S’s defined by (14) and (15), 
namely: 
Q1 


I (¢) — > [210 (¢) Ox (C) = 2m (6) ox, (C)] 


(26) Fea) = De [2mg(C) oap() + 2mplC) of) 


To (Wn Sot) = 3 lig) Eg lttn Syn 6) +g (0) Zp ty, Sy, €)] + Dari my 


(4=1, ..., 013 ¥=1,-.-,%); [20y9(6)], [2on9(C)] and [2,(u,, S,,¢)] being systems 
of primitive periods for Abelian integrals of first, second and third kind respectively 
and m(¢), me (6 ), m, (€) integers. As in the preceding case for ¢ € interior of 
2;,.3(¢) does not intersect any of the branch points 0,(¢) and the number of 
poles (w,, S,) inside 3(¢) is constant. Also there exists a sequence TJ;’(9!) of the 
roots of g(u,¢)=0 such that the roots w,(C), »<7; ($1), lie inside $(€) and the 
remaining roots lie outside. The points ¢; and X,, of (21) are connected by one 
of the algebraic equations q[u(X, ¢),6]=0, p, [wu (X, 6), )=Oord, [4 (X, op =0: 

Since J,==I, we have oe a relation between a sum of transcendental 
functions A(X,_,, X,) and a linear combination of certain integrals of algebraic 
and transcendental functions given by (25). Hence 

Theorem 2.3. Let H= Bit, & &, Xo) be a harmonic vector defined by (1.2) and 
(1.6) with an algebraic associate f(u, S,¢), u,¢ and S being connected by the algebraic 
equation (2). In general the components of Bif, 2 &, Xo) represent a finite number 
of infinitely multiple-valued transcendental functions. Let H\)(X) (k=1, 2, 3; 
y=1, 2,...,n+m)® represent one of these infinitely multiple-valued functions and 
let $1 be the simple closed curve given by (17). Then 


t oj 
De CG. Si Go 1S 4,0) Tl, Sy50) + 
(27) k=1 ve 7;(2) j=1 ve (Soren 
= as Q1 ae 7 ¥ d 
+R — Sf teil) hall) — a0 Rall Ft, 
where A”) (X,_1, X;)= Sai jee me u,(C) are the roots of q(u,C)=0 and C,(X) of 


Oe C= 00 (8) pier (31) are the sequences described above; X, and ¢; are 
connected by algebraic equations. The functions d,(C), Ry, &, and g, are algebraic 
while the remaining functions on the right side of (27) are transcendental. 


8 m is the degree of Q(X, ¢) in ¢. 


30* 
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§ 3. Examples 
1. An Example. Let 
(1) Stull = 4), ee 


For fixed € the Riemann surface 9, (¢) of genus 0 consists of 2 sheets with branch 
points at b,(€)=0 and b,(¢)=1/€; also, b; +b for any ¢ (== ce). The curves Bi (C) 
and $4(C) are defined by 


and 


% 2 4 IL FIC 
S(6) for 7 nO A (EG) for a ha 


2) 


Fig. 1. Projection 3(¢) of curve 3ton the Riemann surface §,(¢) (f=e 


respectively. Let & be the unit circle 
(2) Ve BICC ee f= 2a). 
For €€&, B83 (C) is the line 


ycos#+ zsind =0, %=0 
and 


(3) Bie EX 0). 
The equation of 83 (¢) for €€& is 


ycos?+ zsin? = — sin’, % = cos? 
and 


(4) Bi = E[X|x =cos%, ycos#+ zsind = — sin®, OSS 27'0 
Call B= BB", 

Replacing u by u(C, x), 
() St (ZO oo) CP ee ee 
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In general }i(X) is a 2-sheeted Riemann surface of genus 1 with branch points 
rons Orat 


e,(X) = et, e,(X) ae 
(6) ; 
yy H 4Zz} x —["? : 
o(X = OeNee be Moapone a eed, 


(7?= x*+ y?+-2"). The branch points e;(X) are coincident if and only if 7=0, 
Y4+4Z—0 or +[7?+4Z]}=+r7. The first and last conditions imply X=0 
and Z=0 respectively and the second condition that 2+ y?+2212z=0 and 
y=0. Thus the set G1 is given by 


(7) Gl = E[X |x? + (24+ 1)? =1, y= 0] vE[X|y =z =0). 
Let 31 be the curve 

(8) S1 = E[X|x =coss, y=sins;z=—1, OSs 4a]. 
For €€8 

(9) 3(¢) = E[u| Reu =coss, Imu =cos#sins — sind] 
or for cos #=-0 

(10) (Re u)? 4 (“eee y=. 


31 meets B? in the points (0,1, —1) and (0, —1, —1) and %§ in (1,0, —1) 
and (—1,0, —1) as well as in the two preceding points . Thus X,=(1, 0, —1), 
eNO, =e ve ——(— 1-0, —A), Ag— (0, —1,-—1) and 3" 1s divided ‘into 
four parts with end-points at these points. (Note that X, and X; lie on ©! also.) 
The corresponding points on & are €,=1, €,=(1+7)/23, €,=7, €,=(—1+2)/23, 
C= 1, bg=(—1—1)/24, Cp =—7, Cg= (11/24. 

From (10) we see that the branch point b,(¢) = 0 lies inside § (¢) if — < Te - qt 


4 
b, (C) =cos? —7 sin#, lies inside 3(C) except for 9=0 and a. 


For X €9} 


or % << 2" and outside if BB< 37 or 28 << = On the other hand 


|e, (X)|? = Pars OSS | e.(X)|? = esa 


so that | e,(X)|<14, |e(X)|>4 for —pa<s<sn and |¢,(X)|>1, |e:(X)|<1 
for jcc. Also for sins>0 


4 “A. 4 
cos?s—2 [sin s]? cos s+2sin s 


|45(X)|P = 1+ sin? s 
cos? s+ 2 [sin s]? cos s+2sin s 
fee esl | 


It is found that for 0<s<#a, |e3(X)|<1 and |e,(X)|>1 and for g3a<s<a, 
|es(X)|>4 and |e,(X)|<1. For sins<0, the values of |e,|? and |e,|? are equal 
to the values of |¢,|? and |¢,|* respectively when sins>0. 


2. Example, Continued. Now consider the rational function 


_ p(u,f, S) 
(11) f(u,¢, S) = “q(u, 6) 
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where w,(€) are the zeros of (uw, ¢). If we assume that the surface )? lying in 
E, never intersects $1, then the zeros of Q(¢, X) do not lie on & for X € 3 or the 
zeros of g(u,C) on ¥(¢). Proof: If for some X,€31, there exists a ¢,€ % such 
that Q(f), Xo) =0, then X,€N? by definition of %?. [See (2.8)]. In the second 


= 4 
ae a Te t UB <G< Bue 
4 4 $ (é) 4 4 

S 
9X ege 4 u 
4 4 

0 
t—plane t-plone 


Fig. 3. Transform of 3(¢) 


case X)€ § implies w= (C, Xo) C3 (C) for allé. Ifa zero of g(u, 6) =q[u(C, X),o]= 
C~?Q(¢, X) lies on 3 (€) for some =C,€&, then u (Cy, Xp) ES (Cos ButO¢ (a) 
q(# (Co, Xo), Sol =o PQ (Ey, X). Therefore X G92 which is a contradiction. 
Thus 30 $'=0 implies the conclusion stated above. Consequently in considering 


the integrals and w d ivide 1 i = 
J J i sh e need only subdivide 31 at the points X,(k=1,..., 4) 
BGA At C(j Sade. Ole 
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ht (X) is made into a simply-connected surface by introducing a cut c, (X) 
and a conjugate cut c;(X) as shown in Figure 2. Then ,(X)=n (X)=— 
and m,(X) =0 and since «=1, ge=1 in (2.14) 

Sii(X) = S(X) =420(X)  (@(X) =a,,(X)), 
Si1(X) = S'(X) = + 2y(X) (=m). 


Also if the loops /,(X) around the poles of F(f, S, X) can be drawn so that & 
never intersects them, then by (2.15) 


S124, Sy; X) =+ 2(,,5,,X)  (Q, = 2). 
Consequently (2.16) becomes 


if 


7, a€ 
ee SX) Sse Die A QIG Six) eae (X) ee 
(12) g io 


— 2f(X)@(X) + 27(X)o(X) h=th k=?) 
and by (2.21) 


=> [{d lS, X)e(X) aX +R, (x) dX] — 
— 20 (X) F(X) dX + 2 (X) f(X) aX} 


Since )i,(¢) has genus 0 it can be mapped on to the complex ¢-plane by the 

transformation 
4 2 

(14) t= (#4) Ue") = se ‘ 
where du/di=—2t di€ (#—1)?, S>S()=2t/(2—1), 6, (C) =0>t=0, 0, (C) =1/C +00 
and #=1—-u=co; 3 (¢) > a closed curve 3, (¢) lying on the ¢-plane and symmetric 
with respect to the origin if ¢=0 lies inside the curve. If 6,=0 and b,.=1/¢ lie 
inside (2), then 4, =0 and #,=co lie outside 9, (¢). If 0,=0 lies outside and 
b,=1/€ inside 3(¢), then t=0 lies inside and too outside 9, (¢). In both cases 
t=1 lies inside 9,(¢). See Figure 3. Under transformation (14) f(u,¢, S)—> 
flu), ¢, ¢(@—1)J=fA(4,0), a rational functions of ¢ and ¢ and 


(15) uate ASS c) Bey; ZOIAG 


where >) 7, oe times the sum be the residues of 7, (¢,¢)du/dt in 3,(¢), w is 
the uniformizing variable given by (14). Thus finally 


8 


d da 
(16) I, ={4 is fi u,C, S)d es 3 (0) 4. 
g 3(¢) 45 


Equating (13) and (16) gives the final result. 

3. Singularities of the period functions. According to BERGMAN [2] the period 
functions are singular only along algebraic curves where the roots 7,(X), 72(X), 
7,(X) of the resolvent cubic of the fourth-degree equation (5) coincide. However 
vy; are connected with the branch points e; of the Riemann surface of equation 
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(5) by the relations 


V, = &1 €3 + 0s Ca, Tr, = €y €, —- €5 €3, 1, = €4 Cy + C3 Cy- 
Hence tj = 2 or, (2,4, R=A, 2, 3; 4,4, & distinct) 
imply C= by, CF FG; (2, 7, R=1, 2,3, 457%, 9% distinct). 


These latter relations hold if and only if X € St. Thus the period functions are 
singular along the circle (7). 

The period functions are singular along the intersection of quadric surfaces 
in the following case: 


Theorem 3.1. J/ 


(17) S*= TT [ui —pj,@)) (22, mz2), 
ja 
where p;(C) ts a quadratic function of € (j=1,...,m), then the period functions 


of the Riemann surface R(X) are singular along the intersection of quadric surfaces. 
Proof. Without loss of generality we may consider the case 


ul — (0) =?(Z —a)+ (x —0)C+2*—c 
ul =OZ+4+%0+42*, 
for otherwise we may make a translation of the coordinates x, y, z so that one 


expression reduces to w¢ and the degree of the intersecting surfaces will not 
be changed by this procedure. Branch points occur at the points 


—(¥—b) +13 —«+r3 
2(Z—a) ’ DT aan 


where 7,=(% —6)?— 4(Z—a) (Z*—a) and 7,=x*—4ZZ*. The branch points 
will be coincident for 7,=0 or 7,=0 which lead to circles or in case 


zt 4 
EVs east 


2(Z—a) QZ, 


mare) 


or 


+Zr=ax—bZ+(Z—a) zr. 
Squaring and simplifying we get 
2(Z— a) x a+2b%Z(a—Z)—4Z(Z—a) (aZ* —Zc) =+2(ax—bZ)(Z—a) rH. 


Divide by 2(Z—a), square and divide by 4Z. After some simplification we get 
acx®—abxZ*—bcxZ-+ (B® — 2ac)ZZ* 4+ @Z¥4 2Z2=0, 


which gives in general the intersection of two quadric surfaces as the curves 


along which the period functions are singular. This completes the proof of 
Theorem 3.4. 


4. Singularity curves of higher degree. However, the curves along which 
period functions are singular can be of higher than second degree, as the following 
example shows. Consider the irreducible equation 


(18) S=ul?— 14, 
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which becomes on substituting w(X, ¢) for u 
(19) S?=Z03+ «f?+4+2*0—1 (Zeany 


and which represents a 2-sheeted Riemann surface 3t(X) of genus 1. The period 
functions have singularities where the branch points coincide, that is, when 
the discriminant of equation (19) is zero. This leads after some simplification 
to the equation 


AGZZ* — 47)8 42 (247207 Z*)? —544,22(24— 9ZZ*) 4+-7290Z! = 0 
or upon separating into real and imaginary parts to the equations 
4(3Z2Z* — x*)3 + x (2x? —QZZ*)\? — 54x ReZ? — (2x? — 9ZZ*) + 729 ReZ* =0 


and 
— 2x ImZ? (2x? — 9ZZ*) +27 ImZ* =0. 


After some calculation these two equations are reduced to 


(yy? + 27)3 + x? (2 + 22)? + 2 x (22 — y?) (8 x? | Qy24 92?) DT OP athe) 20 
y z[% (8x2 + Oy? + 927) — 27(z — y)] =0 


so that the curves of intersection take the form 

(aj) 24-4 4222-1 2% (8x7 - 927) — 2727 = 0, =O}, 

(by) yt xty?— 24(827 1-997) —27yY7=0, 2=0, 

RO ys i 28 )2 a (ye tee)" abd (2 Be by) 27 (28 — Oy 22 E ¥t) = 0, 
% (8x7 + 94? + 92?) = — 27(z— 9), 


that is, two fourth degree curves and the intersection of surfaces represented 
by a sixth and a third degree equation. 


§ 4. Generalization of integral theorem to vectors 
whose components satisfy certain linear partial differential equations 


1. Line integrals of solutions of partial differential equation (1). Consider the 
linear partial differential equation 
(1) Ap +F(r?)y=0, 


where A = 02/0 x?-+ 02/0 y2-+ 07/02? and F is an entire function of 7?= x?+ y®+ 2?. 
BERGMAN has shown [4, 5] that there exists an entire function I’(7, t) of 7? for 
[¢| 1: 


(2) I(v,é) =1+ > ?*o™ (9) 
H=1 


such that the operator 
G) yX) =B))=S To.) HX A ae=2 fe, AIX A) ae 


generates solutions of (1) from harmonic functions (X). 
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Let H(X) be a harmonic vector with components H;(X) (j=14, 2, 3) given 
by (4.2) and (4.6) and consider the line integral 


(4) I= (xX) dk=2f [Pe )aeH [Xt —A))-ae. 
aS ~s 0 


If S$! is a curve lying on the sphere x?+ y?+2?=7*, then I'(7, ¢) is independent 
of X and the operator P,[H(X)]| represents a harmonic function L(X); also, 
the order of integration may be changed on the right side of (4) giving 
ae = 1 = = 
(5) I= f L(X)-dX =2f I, dt f H[XU—#)]-dx. 
ay! 0 bays 


aS 
2. Behavior in the large of solutions of partial differential equation (1). In 


order to simplify matters we assume that the Bj-associate of H(X) is a rational 
function of w and ¢: 


_ p(u,S) 
(0) #(,0) = PED, 


p and qg polynomials in w and ¢. Under the transformation (1.3) /(u,¢) goes into 
P(g, X)/Q(¢, X), P and Q polynomials in € and X with Q and g connected as in 


(2.6). Thus 

: Pit, 1—#)X] dt 
H[(1—#) X]= ; 

7) eae J Ol. GA) XTC” 


where & is a closed regular curve in the ¢-plane. 


Following BERGMAN [8] we introduce an R-manifold ©? in £; in the following 
way. For fixed X the equation 


8 B t, rok Se: Q re (1 =e) s 4 SS y(t, Xx) f?— 

( ) ( é ) Gn (t, X) i= gu (t, X) : 

(7, (¢, X) polynomials) defines a Riemann surface T(X) with m sheets = (t, X), 
y=1,...,”, which in general has s distinct branch points ¢") (X), namely, the 
roots of the discriminant 

(9) Dt X) =. 


Equation (9) defines in £; an R-manifold © with s space-sheets ¢ (X), u=4, .. 
The branch line of G? is given by 


So 


oi. 


(10) 6! = EID (i Xx) = 04 D(x Ole 
BERGMAN [5] has shown that if $4 is a path on G— G1, then the genus o=0(X) 


of &(X) remains constant as X moves along 581, 


By the residue theorem 


RPC (ieee An me Rhy: 
0 2 = cr a WY Sh residues at poles inside &, 


g 
where the poles of the integrand are given by the zeros¢=¢,,(t, X) of Q[C, (4—22)X]. 
The residue R(¢,,4, X) at the pole ¢, 1s a rational function of ¢.?and X or an 
algebraic function, sai (t, X), of ¢ and X. Also for fixed ¢ and X a certain 
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subsequence 7(&, ¢, X) of the poles lie inside 2 and the rest outside. Thus 


Pie f= 2x de ) 
: TE Ge ee ee alts 
} , Q Ole) C peT(St,X) ult ), 


where the right side is the sum of certain branches of an algebraic function 
{R,,(t, X)} which is defined on the Riemann surface X(X). For fixed X as we 
have seen in § 2 the interval (0,1) can be subdivided into a finite number of 
pieces, (X) (4X) te (X) nk —ty si PX) such that for ¢e7,(X), a fixed 
set, T,(%, X), of the zeros of Q[C, (1—#) X] lie inside & and also for ¢€ interior 
of t,(X) no ¢,(X) can lie on &. Call 

(12) TG Nas ae tb 2X eae (65 XC) 


uc T,(2, X) 


By the Weierstrass formula for X € 6 — 6}, 


ROG Neg ee Gi ae tO ae "SRO (Cie 
(13) Ne ae 
— 2D [ea(X) Gal 6 X) — feal(X) Gale X)], 


where (t,,,C,,) are the position of the poles of R™ (¢, X) on the Riemann surface 
Z(X), c(X) is the residue of R™ (¢, X) at (é,,¢,) and G,, G,, G are Weierstrass 
integrands of first, second and third kind respectively. If we integrate G,, G, 
and G from ¢,_, to ¢, we obtain 


th a hb z 
Gre X\ db T, (toy ths Xe SOOO 


tha 


i G Le r L,\ 7 
(14) S Galt, X) dt =1q (tea, ty; X) + SEQ (X) 


bea 


th a6 
ij G (t,,.» ch oC; A) dt =I (t,_-1,%; Lose Ce, X) =F Se (Cae » Cee X) 


ta 


where /,, /;, are the values of the integrals of G, and G, respectively over a 
path from ¢,_, to ¢ on the simply connected surface 2’(X), which does not meet 
any of the cuts or conjugate cuts used in constructing &/(X) and J is the value 
of the integral of G over a path from t,_, to & which does not meet any of these 
cuts or the loops introduced around the poles ¢,,; S{)(X), S!"(X), S® (t,, Ce,5 X) 
are linear combinations with integral coefficients of the system of primitive 
periods of the first, second and third kind respectively, entirely analogous to 
(2.14) and (2.15). By means of (14) we obtain an expression for 

tl X 


) 
ROUX) at, 
) 


te_1( 


and summing k=1 to p(X) we get that 


D(X) = -&(X) 
(15) ie GR) Gt 


k=1 tp_4(X) 
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is a linear sum with algebraic coefficients of 

i. the transcendental functions S!*)(X), S®"(X), S (t,,,¢%,; X), which involve 
the period functions of the R-manifold ©, 

ii. integrals of algebraic functions, not involving period functions, and hence 
representing algebraic or logarithmic functions. 

Following BerGMAN [5] we may say that a function which can be represented 
linearly as the sum of functions of types i and ii belongs to class G(S). Conse- 
quently a solution p of (1) which can be represented in the form P,[H(X)| where 
H is a harmonic function with a rational B,-associate is a function of class G(S). 


Wi 


E, (2) 
L 
So 
2, <tst, 
C7 (2) 
- : 
ty = c = t3 
6, (t) 
Surface 3* Subdivision of curve 2 


Fig. 4 


3. Expression for the line integral (4). Let $3! be a closed curve lying on the 
sphere x?-+ y?+ 2272 (0<r<oo) which is defined by (2.17), 9? the surface in 
E; given by 


(10) 9? = B[X 


= # Xs), X(ses 0 sts 1) 
and ¥? that given by (2.8), namely, 


(17) N2= BLX|XE y U xt (¢), 
v=1 CEL 

where 

(18) MN (C) = £[X|u(X, 2) =, (0)], 


u,(¢) being the roots of q(u, ¢)=0. For fixed €, %3(C) is a straight line in Ey. 
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In general 3? and 9? will intersect in segments of a finite number of curves 
si, ..., Gh... The equation of each G may be obtained in terms of the parameter s 
by substituting (2.17) for X into u(X,¢)=w,(¢), separating the latter equation 
into real and imaginary parts and eliminating ¢, thereby obtaining a relation 
between ¢ and s. Solving for ¢ in terms of s and substituting back into (2.17) 
the desired equation is obtained for GF. In general G will be defined for s in 
a finite number of closed intervals of s and possibly at a set (assumed finite) 
of isolated points. Let 2’(s) be the line on 3? defined by X=#X(s), OS#X1. 
The curve G} may coincide with some @’(sy) for ¢ in an interval, in which case 
Sp Will be one of the points of subdivision of s. Otherwise suppose that G} inter- 
sects each Q’(s) a finite number of times. Since G; is a curve and Q’(s) varies 
continuously with s, the number of points of intersection of G; and Q’(s) will 
be constant for s in non-overlapping intervals with jumps at the end-points 
of the intervals. Thus s may be subdivided into intervals in which G} intersects 
each 2’(s) a constant number of times. Repeating this procedure for ©} (k + 7) 
we obtain a subdivision of s: 0<s,<---<s, such that for s€S,=(s,_1, $,) 
there are p, points of intersection of {Cj} and &’(s). Let Y be that part of 1 
for which s€ S,. Then for X € 3% (15) becomes 

Px te(X) 


(19) es R(t, X) dts 


Rk=1 ty 4(X) 


ss 
Consequently the line integral fy (X) - dX is broken up into: 


2 aa me = 
(20) Lye aA, 
#=1 x=1 wETys(2) Bz 
where py” is a vector whose components are functions of class G(G) and 
(X) 
(21) yl) (X) = i I'v, t) R,(, X).dt. 
th 


4, Relation between the line integral (4) and certain functions of one complex 
variable. Substituting (7) into (5) and interchanging the order of integration 
in the last two integrals, we get 


1% 
fees dt ¢ plu—e)u,] 
vaafreyat {4 | q[(1—#) w, €] 
0 2 (0) 


1 
2) I'(7, t) ac Pv, ¢) 
2 f 1? Ki i! ¢ | q(v, ¢) ay; 
) g V5) 
where 9 (¢) is given by (2.23) and 


(23) VC, 4) =£[vlu=u(X,¢) (1—#), XEO'. 


We assume that 9’(f,¢) is a closed curve except possibly for a finite number 
of values of €. By the residue theorem 


»¢ 
(24) f Ses -40=2r 0,01, 
S(6, 2) 


(22) 
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where 7, is 2% times the residue of f(v,¢) at the poles v,(¢) which lie inside 
S(¢, t). As in § 2 for each ¢ we can subdivide & into a finite number of pieces 
@,(t) with end-points ¢;_(¢), ¢,(t) in each of which the set of zeros of q(v,¢) 1s 
constant. Thus 


(25) 2 f AM) ae 


In the considerations of Paragraph 3 we could have subdivided (OS#<1) intoa 
finite number of pieces 1, = (t,_1,¢,) and 3’ (¢) into a finite number of pieces = 
(X,_,(#), X;(6)) such that for t€7, and X€¥;(¢) the number of points of inter- 
section of $2 and the curves {Gj} is a constant (3’(¢) is the intersection of 3? with 
the plane ¢=?). For each ¢€t,, X;(¢)€ 3? implies that there exists a y and a 
C(t)€& such that X;(é) EM (C(é)). We assume that X=X;(f) intersects each 
M1 (¢) either in a finite number of points or coincides with 9} (¢) for ¢ in a finite 
number of subintervals of t,, which case gives one of the points of subdivision 
of 2. In the former case since X;(¢) moves continuously with ¢ and Ni (¢) is a 
straight line, t, can be subdivided further into a finite number of non-overlapping 
subintervals on each of which the number of points of intersection of {Cj} and 
93 (¢) is constant with a jump at the end-points of the t-intervals. Reletter the 
subdivisions into a simple sequence Tt, (\=1,...,¢) such that for f€t,, g, of the 
points X,(z) lie on 9 (C), say €, (4), ..., €4, (¢). These points also lie on & and hence 
subdivide it into qg, pieces &; with end-points ¢,;_,(¢), €;(¢). If the v,(¢) are the 
roots of g(v,¢)=0, as in the case in Paragraph 3 there is a sequence T,;(9) of 
the v,(¢) which le inside the domain bounded by Q’(t,¢) while the remaining 
v,(C) lie outside. Thus (25) becomes 


(26) 2 yy faa [rb.0.4. 


Equating J in (20) and (26) gives us the following 


Theorem 4.1. Let p be a vector whose components are solutions of (1) which 
can be represented in the form P,[H(X)| where H is a harmonic function with a 
rational Bs-associate plg. Let S% be the R-manifold on which is single-valued. 
Let .* be a simple closed curve situated on the sphere E[x®+y2+2=72, 0< <ool, 
satisfying conditions noted above. Then 


Re Pare 
(27) pia De ae eS 
m=1 R=1 wETyp (2) Xp1 


1s equal to (26), where v,(¢) are the roots of q(v,¢)=0 and Tp (X) and T(S) are 
the sequences described above. 


Presented to the International Congress of Mathematicians, Edinburgh, August 
1958 and by title to the American Mathematical Society, Eugene, Oregon, June 
1959. 
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Ouadratic transformations of generalized Legendre’ s 


associated functions for special values of the parameters 


B. MEULENBELD 


Communicated by R. TIMMAN 


1. Introduction 


In previous papers [1], [2] Kuipers and the author gave two linearly inde- 
pendent solutions P””(z) and Q;”"(z) of the differential equation: 


aw dw m* n* pe 
(1) (1 2") Ag 2% dz {Re (# 1) 2(1—2z) sate Oe 


specified either in the form of integrals over double circuits, or in terms of hyper- 
geometric functions valid in restricted areas of the z-plane in which a cross-cut 
exists along the real axis from 1 to — oo. In other papers [3] and [4] we derived 
expressions for the particular solutions of (1) 


Fr) Pet PG Pa ea oe 
OF @, OA"), OMG, Brl—H, OFm(—2) 


in terms of Py””(z) and Q3"(z). 

By means of the well known linear transformation formulas of the hyper- 
geometric functions the Py””"(z) and Q’"(z) are expressible in several ways in 
the forms 

Be (2) AvP tay bye Gh apy PA Ga oa eee ee ad 
and 
grt (2) == Ag FP {ay 033 Cg; Ch AVP tag, gs ge Che AG 


where ¢ is a linear function of z which depends on the choice of the transformation. 
In [5] we carried out these expansions. 


If the a, 6 and c in & {a, OS Os zh are unrestricted, there exists no quadratic 
transformation of the hypergeometric function*. If and only if + (14—c), +(a—b), 
-+(a+6—c) have the property that one of them equals $ or that two of them 
are equal, there exists such a quadratic transformation, this being given by the 
formulas of GAuss and Kummer. For unrestricted values of k, m and m in Py””(z) 
and Q;”"(z) the conditions for quadratic transformation are not fulfilled, and 
therefore it is not possible to derive quadratic expansions of these functions in 


* I.e., a transformation of F {a,b;c;z} in functions F {a’, b’; c’;y}, in which 
the variables z and 7 are connected by a quadratic relation. 
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terms of single hypergeometric functions. We have deduced expressions for 
Ber" (z) and Qf’"(z) in terms of series of generalized hypergeometric functions 
or MacRosert’s E-functions. These results will be published elsewhere. 
However, when the parameters m and u satisfy special conditions, quadratic 
expansions can be derived for the Pf"(z) and Qz°”(z). In the present paper 
I shall deduce expressions for these functions as finite sums of hypergeometric 


functions with quadratic forms of z as variable under the rather strong condition: 
n—mM 1S an even integer or 0. 


In §2 these expressions for B””(z) are deduced. In §3 I shall prove a 
relation between Q}”(z), PB,” ~”(—2z) and P,””(z), by means of which the ex- 
pansions for Q2””(z) are deduced in § 4. 


For m=n these expressions lead to the well known quadratic expansions 
for the associated Legendre functions (see [6], Chapter III). 


2. Expansions for P’’"(z) 
The starting-point is the relation {/, 11]: 


mn 1 Z=-4 r; n—m n—m . he les 2 
Ne Vere Tam ; wet { ae +4, k+- ae {—m;- \ 
ae 


valid for |z =4| | < 2. On account of the relation 


2 Fiat pa ee TO 


Spe cle 6 ey 
T@reretn n,b+n;c+n; 2} 


when n — m is an even integer or 0 we have 


in eta i LED 
Pee eNom (5-4)? rk 2 +1) ae 
h ea = eee 
(z—1) 2 r(k+ Z * 44) r[1— 2") 
3) n—m : 
5 d ae Fk ees BoA wen “Sh. 
ig - 
Since 
m : z+1)2 ie 
(4) 2) T(—m) (4) FUR +4, —k-1—m; ei 


where P?”(z) denotes the associate Legendre function of the first kind, the 
expression (3) may be written: 


ea 


(z—1) 2 r(k+ So +1] 
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2.1. For |z| <1 we have 


als Cs 
2m ni? (z2@—1) 2 
Pi” (z — — 2 es . 2 x 
Pt) Spe Pe eae ee 
|e ee MR css DESMA, TA 2} 
et Be ire C m 
eS 5 5 Maa ; rises ac 


T@h—tm+4)T(—4h-¥m) 
fear he wal eA roe )) 
cana ai Ws ea ap as tees f 
Hence (5) can be transformed into 
& Fil ee ee 44) 2" mt 
Pp” n (z) (a4) 2 ‘ = = = xX 
We n—m 1 m+n 1 1 mtn 
cals P(e D 1) r( ae 4 isin 4 +1) 
n—m ne 
a -4 1 mtn 1 min,1 1 a} o 
ns ma | +1) PY 2 4 Nae 4 +4; 5; 4 
6) Wa 
te Lew ; nm+1 2 
(Zar Aye r(k Z +1)2 a : : 
Le n—m 1 m+n 1 es fl m+n \ 
Raa pee yee 4 apie a eae 
nam mare 
ay eer ar 1 m+n 1 1 m+n Nice <r] 
x aoe e+ Peed ri are ree its all. 
2 


We find, on carrying out the differentiation and employing the relations 


ORIG Glee ee 


’ 


sin 7 Z 
d 
an T(22) = 2! 49-31 (2) (2 +4), 
that 
(7) Bm" (2) = A,(m, n, k, 2) + A,(m, n, R, 2), 
where 


(—1) 2 9-71 r(* = +1) r(a— ">" 44) (z+1) 2 (2-1) 72 
= r(@E")r(s Bar 1) (4 Es Hi) P(—a a : 
ese Yr Y 
(8) x 2 Z rs or) le, 
: el ay-*er( sh a 4y a\0(5r arn yas Is 
7-0 P—2g-+4) Pgs) F(5 +74) 
xF Sh me ty q 5k mem ty qtsigtr qh, 
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P(A) r(a me 


2 


; 
is : 5 P(n—r 1) e+ A)r( aye {2( 2S 1] eee } 


; sl 1 Mn , 
% 4 Ly a+3)r(Sh a LY oe 
ms Py—2g-+1)Mg+1) P(5 +1—4] 
1 1 1 
x F{ eae matty Rares Rory mt ty g+4;3+7 gah; 
these results hold when |z| <1. 


men 
2.2. If we substitute in (5) for B, 7 (z) the expression 


ol is mn . 1 m+n , \ 2 
4g Ea 2 ret. ype \, 
valid for \4—2?| <1, #z>0, and follow the same argument, we obtain: 
Peet (2) = Bim, n, ; 2) 


MS $1) amet) ® (+4) 2 
Z x 
m+n a 1 m+n ft 1 m+n 
ra") ri 2 har(ae 4 +a)Fl ies 
a 
I'(n—yr) (2+1)’ 
“ » im—m os 
7=0 P| 2 =) 
(41) 
[z] (—1)! (2224 (52 dD LEST q+5)r( Lay a) 
x B 4 2 2 4 y. 
= P(g+1) Py—2q41) F[1— ""™ +r—4) 
1 mtn 1 1 m+n ‘ 
ame 7 f Slgaa Bk 5 +r —q; 
nl mn , 
which is valid in the same z-region. 


ee gq; 4 2}, 
2.3. Applying to the hypergeometric functions in (11) the quadratic trans- 
formation 
ee LC el (O— a) = { et | 
f i x, b 
ee (6) Fe=a) (—=2)n °F a,4— c= a3 4 ba; — 
Pc) P(a—b) (_ »-» J ! 
E F(a) Flea) 2dr OMA re, 


\ » I 
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b;1—a+b,; At |arg(—2)| <2, 


31a 


x 
r(h- ae +1) r°( h =| 
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we find after some calculation that 


(13) Pm" (2) = C(m, n, k;-2) + C(m, w, —k —4, 2), 


(NEE it eee +4) r(;— 7S Hi) r(—2- 3) x 
a aa Pe 1 m+n 1 
oa arene eens 
Jay n= 4 7 ; or $ a 
ey ees : 3 Nhatrgea eae 
x Th mye on rs +4) 
: 1 + : 
ie (= 1)? (22)? 4 (2+ 4)" (2-1) 144 P(e min Bid >] 
P @+1) P—2941) 
d= 
1 m+n | 1 tg SE lig pile : \, 
XFL & zl eae ee 2? rai 4 rede 2°? 4—22 ? 


by analytic continuation, this formula holds when |1—z?|>1. 


2.4. If we transform the hypergeometric functions in (11) by means of 


(15) F {a, 636; 2}=(1— 2)" F fb, ¢ — a3 6; 4 \, 


—— 
a 


we obtain 


ee (A\ ee (He 00s 2) 


(ay (a ™ +1)r(k or, 7 (1) (ep 1yre Ata 
x 
m+n n—™Mm 1 m+n 1 1 m+n 
r( rt) r(r 5 Hr (Sk ; ba) rl k 
gee we 
(16) . Ny a oe ae 
v=0 r| —r+4] 
To i A 
Sine idee =" sides ks ral xt ae ee 7 
x 
= Pg+t) P—2g+1) P(t =" 47-4) 
FE " shes m—-_n ee ul m+n F m+n : 1 
x { 2 a 4 T RP 2 k 4 Yr qg; 1 2 med Gai 72 ? 
which holds when ]1— 5 <A, Jus O- 
z 


2.5. On applying (15) to the hypergeometric functions in (13), we have 
(17) Bye" (2) = G(m, n, k, 2) + G(m, n, —k —A, ay, 


where 
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m n—-m =, — 
Citra tart 1 (75™ 41) r(r— 5™ 41)x 


ees ("SE") r(e+ ae is +1) P(2 peuple ae a|x 
2 \ 2 2 ga.) 
1 n Oy ee ath ema 
x r(—h- 5 (=f2Getyre sas Py) ‘ 
mtn 
x P(—k— : = joo T(7S™ 141) 
(18) . ; one , 
Ce) (<1) (2ey'-24 (2-41)! (eater (5 ie fey ot 
% 2, T@+1) 0 72941) i 
q= 
1 mtn 4 1 m+n d Shee hi! 
xF {ER Eh Sh EM tg ttskt os ah, 


which holds when |z|>1. 


2.6. Furthermore on applying (12) also to the hypergeometric functions in 
(16) we obtain 
(19) PY” (2) = H,(m, n, Rk, 2) + H,(m, n, R, 2), 
where 


H,(m, n, k, 2) 


ae jFni(dat = a) on (=™ ™ +1)r(r— "5" 44) 
: x 
m+n n—mM 1 ee 1 : ss 
Behe tte ge ealtlag 
eae er ee 2p 2om 


Cy as (et ya ee x 


re nae: [2] (ayer page 
: 2 es ne 7 Di  apatlanaes 


2 4 oy 2 4 2 2 fp 
and 
a ee he ,_ mam 
Ate bem a ee ee +1) 
Page m+n nn. m+n ee 
r( : )r(r +1) r(k— +1) 7 ph 
Geetha CER) waa te x 
eS [z] 
I'(n—r) (yh (a7y 74 Gat)! (21) x 
(21) 4 > T(q+1)x 


1 m+n 1 Aer ents isn 
xI(>h 7 +y ats) P( me 4 Ee a) 
xP(r—29+1) (5 +79] 


Z2—1 


2 
x F{-fk- i +74, hp et Stra if 
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the upper or lower sign of the exponential being taken according as §z>0 


ge 
a |< 


or <0. Formula (49) is valid when 


3. Relation between Pj,” —™"(—2z) , PR’ ™(z) and Q7"(z) 
Theorem. When n—m is an even integer or 0, we have the relation 


L mim 


—nmim m,n — __ Whe —— 
f Oe) 2sin wk a 
(22) 


Pa aol z) Pas ermtk Peony) F 


a tiga a itt Behe See 
WS (2) Qn r(et ™" 44) r (a or +4) x 
Y {gr (2 — 2 sinmme*i™ on (2) [4, (48)], 
4 ma n—m 
i 
ey 


n—m 


x{ = sin (A mn) we tm OM (2) | erui(i+ “3") Pee (2) [4, (21) ] 


and 
q n—m 
orm (— 3X ett (+ k—m-+n) D (e+ > +1] 2 
k z) an—m r(h n—m Rk (z) [4, (25)], 
+41 
we find 
r(a— m+n +1) 
Bot (8 oy ap a -|2 SU OR ater sin (4 Tite a 
x 
rat si 7 mu ( ) | 2 


pesie red. een Eres : 
<cage +en** sin nae} + eF ** Pym (2) 


P(e mm +1] 
= — te 
r(a+ a +1] 


Bp fal, - 4 . . 
\ —e mim Oe" (2) sin whet tim4 oF ih pymn ()h 


By solving for Qz””"(z) we have (22). 
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4. Expansions for O}" (z) 
It has been shown in [5, (8)] that 


ror" 


on—m— ur (h+ 3 aaa +1) (k— is” +1) T(—m) ; 
= ers Fa (a4) 2) (21) 
r(h+™ 5 )r(e- : +1) 
XF {—k—2S" = 4; 14+m; I 
(23) 1 n m 
+ = Lim) (g=-1)2 (¢-—4) 2 x 
x FI k- aS h “= Ae Aa is SI}. 


Employing the well known relation 


as ies eee ont ma ecre gy 


ene 21" (4 — zt? ™ Bla —n,b —n;c —n; 2} 


and (4), we find that 


F{ Fee 4 | 


2 2 : eae | 
(24) n—m (e+1)" Ean ca mee _ mtn 
= py | : ¥ ee P. 2 
2 I I m) (z— yea eee ( z+! k a) ? 


and, on account of (2) and (4), 


Fi iB poem RY n—m | m: tz | 


2 z= 
- LS (ae oh 
From (23), (24) and (25) we obtain 
preter es 
(26) “| ‘ — ao (Nee ie a+ 
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4.1. Since —m=—n, we may apply (22) to (7), and we find that 


nb E ey 


mexknmim 


A,(—n, —m,k, —2z)+ 


enim mm) = — SE 
2sin awk (h-™E" +1) 


| A,(—n, —m, k, — 2)} eFih fA (m, n, k,z) +A,(m,n, k, 2}] 


where A,(m,n,k,z) and A,(m,n,k, 2) are defined in (8) and (9) respectively. 
From (6), (8) and (9) it follows that 
A,(m,n, k, 2) 


n—m no 
ed ae r(k— +1) 2 IU 
ti mn n—m 1 mtn 1 1 m+-n 
ot Hee ae oe 
aa a eet { mtn 1 mtn 1 1 2 
2 AH Bic! 
«ee [+4 F{-ta—-7E", Se- "4S SA 
A,(m, n, k, 2) 
n % We n -- 
(2+4)2 ri +4) ae ; ree. 
au n—m 1 m+n 1 1 m+n 
cot Te = Hr er oe 
nu—™m 
a Bite 1 mtn 1 1 mtn 3 
2 pee b Bs asl eae eae 2 


By changing (m,n, z) into (—n, —m, —z) we find after some reduction: 
A,(—n, —m, k, — 2) 


1 1 f : 
robe MM A) (La mE a 
= LS EAS : —A,(m, n, k, 2), 


A,(m, n, k, 2), 


mime | GE 
rte Et 


—— A,(— n, — m,k, — z) — e**** A, (m, n, k, 2) 


sin (} k— are ao 5) is 
a ern A,(m, n, k, 2) 


- x) 


i Ae L 
gre a 19 2 
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In a similar way we find 


P(e+ +4] 


P(r oe +4] 


A,(—n, —m, k, — 2) — e*™** A,(m, n, k, 2) 


e 2 ‘ 
= 1 ean sin wk A,(m, n, k, 2). 
sin (> k+ ate 
so that 
ere Or Ke) It sl eee era as a A h a: 
é a ee : 1(™, n, k, 2 
cos & k+ meen 
(27) ateinen fay 
oA e 4 2 
Nl aioe Ae A,(m, nt, R, 2). 
sin (| k+ 4 4 Jn 


4.2. To derive an expansion of Q}””"(z) in hypergeometric functions in (1 — 2?) 
the method used in § 4.1. fails, the formula (10) not being valid for R(z) <0. 
Therefore in this case we shall use the more direct method. 


From (26) and (10) we deduce 


ett OM (2) 


pea n areal m+n 
Zils eRe eee BGM igs 1) 
; n—m- AL me =, mtn m+n 
r(h+ 5 +1)sin ma (z=1)? r(k— me tart 5 ] 
ee i es 
as oa 2 4 Be 
Ob 
1 m+n m+n AN 
Res ae ets 2 at 
n—1 ae eas 1 1 
| 2 m+n R222 {(e ! 1) : ms : a r Be 
r| 2 Oe aes 
1 h m+n mtn Rei a) 
2 7 ae 2 
On carrying out the differentiations we find that 
—nmint —Qm,n a Te 
8) AG oe 2sinm x 
2 mtn 
r(k+ mon +1) $ni ee 
x1 B(m, n, k, 2) = — B(—n, —m,k, — 2)>, 
pe) 
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where B(m, n, k, 2) is defined in (11). Formula (28) holds good when [J1—2#4| <4 
and {tz> 0. 


4.3. On substituting (13) in (22) and following the argument used in § 4.1., 
we find that 


— in 


2 Se 
~ia) a 


(29) enim Om (a) = EU 


Mm; n, Ry 2), 
sin (a+ 


where C(m,,k, 2) is defined in (14). Formula (29) is valid when |1—2z?|>1. 


4.4. On applying the transformation (15) to the hypergeometric functions 
n (28), we find that 


—nim ~—m,n oa aU 
4 Qk i ~ 2sinm x ‘ 
O : 
(30) rales m+n +i)e $nih 
x 4D(m, n, k, z)+— 2 —— D(—n, —m,—k—1, —2)}, 
ren EE a) 
where D(m, n, k, 2) is defined in (16), and where |1— e <1; 


4.5. Substituting (17) in (22) and following the argument of § 4.1., we obtain 


—nim -—m,n ies ene eee ; 
61) CO ta ae 


where G(m, , k, 2) is defined in (18), and where |z|>1. 


4.6. In a similar manner we find, substituting (19) in (22), after some cal- 
culation 


, cm 
enim tM (2) i oat glee a AOE (m n.R z) 
ms mag 9 ’ ’ ? 
sin (5 a) 
(32) , 
RL ei gas: a 6 a+ 4 
+ + ——__________ H,,(m, n, k, 2), 
A cos (| k+ a 
2 4 


where Hy (m,n, k, 2) and H,(m,n, k, z) are defined in (20) and (21), respectively, 
and where |- =| ai be 
aaa 


[1] 


[2] 


[3] 


[4] 
[5] 


[6] 
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On the Convergence 
of the Rayleigh Quotient Iteration for the Computation 
of the Characteristic Roots and Vectors. V 


(Usual Rayleigh Quotient for Non-Hermitian Matrices 
and Linear Elementary Diuvisors ) 


A. M. OSTROWSKI 


58. In the case of a non-symmetric matrix, use of the generalized Rayleigh 
quotient appears to be particularly conformable to the problem, and in any 
case much more so than the use of the Rayleigh quotient in its original form. 
However, in some computations carried out in the mathematical laboratory of 
the Ramo-Wooldridge Corporation a routine implying the usual Rayleigh quotient 
has been used by Dr. FRANK and turned out to converge quite well. It appeared 
of interest therefore to discuss the application of the usual Rayleigh quotient 
to the determination of the eigenvalues of an arbitrary (not necessarily Hermitian) 
matrix, and that will be done in this part of our series of papers. The main 
result is that the convergence is indeed essentially quadratic, if to the eigenvalue 
in question correspond only linear elementary divisors. Thus the result appears 
to be weaker than that achieved in the case of the generalized Rayleigh quotient 
iteration, where we get a cubic convergence. However, if we denote in this con- 
nection the computation implied in the solution of one linear algebraic system 
of order ” as a “‘Horner’’, the application of the generalized Rayleigh quotient 
implies spending two Horners, while in using the ordinary Rayleigh quotient 
iteration we can achieve with two Horners two consecutive steps and therefore 
get the fourth degree of the error. From this point of view, even in the case 
of a non-Hermitian matrix, the use of the ordinary Rayleigh quotient iteration 
appears to be not only permissible but even advisable. 


59. For an m Xn matrix A, we consider the iteration procedure starting with 
a column-vector &) and an approximation A, for the eigenvalue and continuing 
according to the rule 


(212) Gp = (AA ele rene ie ee 
* — ? ) 
(213) 1, = oe | 


This paper was sponsored by the Office of Naval Research and was prepared 
in part under a National Bureau of Standards contract with The American University, 
Washington, D.C., and in part at Stanford University, California. I am indebted 
to Mr. G. CuLLER and E. V. Haynswortu for discussions. 
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Let o be an eigenvalue of A of multiplicity m, to which correspond only linear 
elementary divisors. Then we shall prove 


Theorem 1. There exists a positive 6=6(A), depending only on A, and a 
matrix H, depending only on A and o, such that if we have 


(214) [Ay — o| S 6, 
and if for an eigenvector 7 of A from the right belonging to o the relation 
(215) [50— | S 6 || 
holds, then the sequence i, tends to o. 
Theorem 2. If the sequence A, tends to o without any of the i, becoming equal 
to one of other eigenvalues of A, and if & satisfies (215), then we have 
(216) A414 — 6 = 0(A, — 0) (v > oo); 
and even, for any fixed integer p, 
(217) Ay41— 6 =0((A, — 0) (Agee ONece (Ape 0) (v + 00). 


Theorem 3. Suppose that under the conditions of Theorem 2 all eigenvalues 
of A are real and correspond to linear elementary divisors. Suppose further that 
the following condition is satisfied : 

(a) o ts not an arithmetic mean of two different eigenvalues of A. 


Then we have either for a certain constant IK 


(218) jus — 9 =0(K TT (4,—0)*) >) 
x=0 

OY 

(218°) Goa —>g+0 (v —> 00). 


If on the other hand the condition (a) ts not satisfied, either (218) holds or the 
relation 


foxe} A, heey oO , 
ee) aaa 


either for even v or for odd v. 


60. Since, as we assumed, to o correspond only linear elementary divisors 
of A, there exists a non-singular matrix S such that we have 


(219) SASt=clh, t+ WV +cl,_,), 


where V is an (n—m) X(n—m) matrix in the Jordan canonical form, all eigen- 
values of which are +0. 


Put 

(220) Oia Gt Ait AS OFA ped) 

and 

(221) ead [6s (Opes ee Peas 
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If we multiply (212) from the left by SP,_,, we can rewrite it, using (219) and 
(221), in the form 


[Bale HV Otte BAS $, = 9,4 FS p13 
or, replacing » by »+1 and putting 
(222) P_,Sé=¢, | (v=0,1,...), 
(223) [OL (Ver Os] gm) eet Open | Males) 
Express ¢, according to the decomposition (219) as a direct sum 
(224) Cape, 


where f, is an m-dimensional and «, an (m —m)-dimensional vector. 


Introducing this into (223), we obtain two equations 


(225) OL owe = OB 
and 
(226) (V+ 6, Ye Ay 47 ae 0, Hy 


But (225) gives B,,,= 6, and therefore we have f,=f (y=0,1,...), putting 
Bo=f. By (222), the decomposition (224) now becomes 


(227) Cy = Py yl Soe Boley. 
61. From (227) we have 


P,_,P_ é& &, = (6* + of) St48 UB -ho,) 
or, putting 


(228) De Oe 
where D is a positive Hermitian matrix, 
(229) [Feo ee—— Ge Dee. 
Further, by (227) and (219) we have 
JR sr 4 é,—o|PiP a =P, P17 (A —ol,) &, 
=o, 8 1 SAS(A— ol) ot, = DIO, V) Gs 
where O,, is the (m xm)-matrix consisting of zeros. We have finally 


(230) |\A-a 


(EX AE, —o EF E) =C¥ DO, + VNC. 


v 


If we subtract o from both sides of (213) and multiply the numerator and 
the denominator in the difference on the right by | P,_,|?, we obtain the expression 
(230) as the numerator and the expression (229) as the denominator, so that 


231 Ale wap ct DO». + V) oy 
(231) y= ae 
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62. We now write D in block form corresponding to the decomposition (219) 


d 


(232) b=(5" ae 


21 Do» 
and then obtain by (227) for the denominator in (231) 


ct DE = (pt ton) (5 5) B+a) 


: / 


= (B* Dy, + oF Dey + B* Dy. 4+ 05" Dy) (B + «,) 


= PD p+ oo D1 8 + B* Dy a, + oF Doo %,. 
Putting 


(233) |B)? Do =B* D1 8, |B)? Dy = a," Do B+ B* Dy 2%, |B|? De = a" D320, 
we see that the denominator in (231) is 


[Bl? (Do + D, + Di). 


As to the numerator in (231), by (232) we obtain 


i W\ = OF D5) 
D(On ) Ta Ko Dae v | 
and 

. OD Z 
* * 1} 
B 0) (5 D5 V 


Therefore, if we put 


(B = a, ) = B* D,,V a, + a, Doo Ue Bs 


22 


(234) |B|?N, = B* D,.Va,, [plzNG a Dg2V «,, 
we have 
eek N,+N, 
(235) 0 Delp een: : 
On the other hand, putting 
(236) Viale 


we have from (226) 
VND eine is by W) Oy = oy Ay 
and replacing vy by y—1, 
(237) Oe Op4 Meare Opt ality 16 
63. From (227) we have 


56, = 6%: 


Denote the square roots of the smallest and greatest eigenvalues of the positive 
Hermitian matrix S*S by A, and A,. Then putting S7=—t, from (215) we have 


Ay|o— | 2 |So— 7 


Ay 
(238) ko— t| 5 41 5 |x]. 


0 


“) 


[7] = All, 


476 A.M. OsTROWSKI: 
Putting »—m=l, let the vectors a and f’ be in their corresponding spaces 
B=) (Hrs vee adel Oy HA Va eee, Ws): 


and for the m-dimensional eigenvector t of (238), 
TSM orate Overs Olle 
Writing 
(239) er 
0 
we have then by (238) 
l m 
DAD tees 
A=1 “=1 


4, — he e2|r|2, 


1 m m m 
ou (Male Dy lel? Dial? ada eae) te EL 
A=1 ri “u=1 = 


ax 


> 


A 
x 


Hae al oe Pee Kj Size yp Bey) 


x 


further, if we apply the Cauchy-Schwarz inequality and assume e<1, 


? 


p? + (1 — e*)|r/?S 2/8] |e 
|B] 1 p _ Vie [rt] 
25 pet 


pvi—e — 2\Vi—#] 2] * 
|B |? = (1 — &) p? = (1 — e*) (|ao|? + |B?), 
(240) eerie 
64. Now we put 
_ || 
eee ae ai 


Since D in (228) is positive and therefore also D,, and Dy, in (232) are positive, 
from (233) and (234) we have 


DZzby, DStiy,, D>, 
IN| Spey, No] S Psy, 
where 9, 21, 2, fz are positive constants depending only on A. Thence by (235) 


(242) |4,—o| =|6,|< BatPsr y, 
Po—Pi VY» 
provided y,< Po. 
Pr 
Put 
(243) 4) = Max (1, 5 (bs +hs)). 


Denote by c the exact bound, ||W ||, of the matrix W in (237), t.e., the square 
root of the maximal eigenvalue of WW*, and put 


(244) Cs = Max(t 20)5 
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Finally, put 


= 4p 1) 
(245) y= Min (3 2 B,? al 


We denote by x a fixed number such that 0<%<1 and assume that for a y>1 
we have 


x 1 
(246) Se peepee 
Then it follows that 
WU + 8, . WAS E18, [NI 20, 
u=0 


and therefore from (237) and (241) 


lieg| Seem Mie = 


(247) Vy Sy 6, 1 VeSSs ¥Y Aas 
Further, if we assume for a y=1 that 

(248) Vr So, 

it follows from (242) that 

(249) = pe + Ps) ¥y Sa Yp- 


65. We now make the assumptions 


(250) Mo = Cy, 

(251) Jools 

hence we shall prove that 

(252) a gee ae (v==41, 2d; 
(253) [6,327 eee 


2 


When v=1, (252) follows from (247), since 
M1 SX*Yo SX. 
If we now use (249) when »=1, we have from (248) 


| d,| Sqn SqqxS— 
Co 


> 


since by (245) ¢,¢3S —. Assume now that we have 
2 


(252°) EE is, 
and 
1 9v—2 
(253°) |O,-a| S 
2 


Then, since (246) is satisfied, we have 
Qqv-2 


is Ve a 
Vy» = 1) é C3% Saye 


2 
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from (247), and (252) is proved. Now since (248) holds, we can use (249); com- 
bining it with (252), we obtain 
\9,,|Sa7S es ee 

and (253) is proved. 

We wish now to choose 6 in such a way that the conditions (250), (251) 
are satisfied. But from (240) and (239), using (241) for y=0, we have 

Ay 6 

oe JAlys 


yji-e = Sy ee A 
1A aw 


Vos 


The right-hand bound is certainly less than cs if 


2 


Ay so) =< 2 
i 6 <= 


ct, OO< 2 Ao 
0 3 a eet ee 


1 

since by (245) cs > On the other hand, (251) follows from (214) if OE 
2 

with 0<x<1. Therefore it is sufficient to choose 


60) = Min ( WB a) 
Cc 


pee Sar gm 
yp 3 eA 


to insure the formulae (252) and (253). Theorem 1 is now proved. 


66. To prove Theorem 2, assume first that the inequalities (250) and (251) 
hold. Then we have from (249) and (247) 


| 0, | = ay, Cy CoVy—-1 | J, 3| , 
and this proves (216). Using (247) # times, we obtain 


Operon io, A, 1= 0; 105 2 


n2 A 
OS — S¢,| 0, 10, y EOS ay yee 


and this proves (217) since y,_,—0. 


If the initial conditions (250), (251) are not satisfied but 6,0, it is sufficient 
to prove that for a certain index k we have 


|| S—, lvel Ses. 


Since 6,—+0, the inequality (246) with x<4% certainly holds from a certain » 
onward, and since from there on (247) is valid, y,+0. Theorem 2 is completely 
proved. 


67. In proving Theorem 3, by what has been said in the proof of Theorem 2, 
we can assume without loss of generality that the conditions (250), (251) are 
satisfied and therefore also the relations (247) and (249), since otherwise we 
could replace » by »--k. 


By (241) and (233) we obviously have 
D, =Oly,), D, =O(y;), 
and therefore in virtue of (232) Dy +D,+D,—+D,. Thus by (235) ‘since: D5 0; 
(254) D(A, — 0) = (1+ ¢,)(N,+N,), 2,0. 
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From (234), N,=O(y%), and since by (247) 


(255) Yr S|P,_1| 70. 
finally from (249) 
(256) 0, = 0 (63" B_4), 
(257) NOt” P24). 
We wish now to discuss N,, given by (234). Introducing the vector y by 
(258) ae Sar B* DioV, 


we have for the (7— m)-dimensional vectors y* and «,: 

(259) iN SEM as 

If the vector y defined by (258) vanishes, then by (255) we have from (254) 
i, — a =O(N,) =0(c3" P?,), 

and we have the assertion (218) of Theorem 3. Discarding this case, we can 

assume y* + 0. 


68. Under the assumptions of Theorem 3, the matrix V in (219) is a diagonal 
matrix with the elements o,—o along the diagonal, w=1, 2,...,%—m, where 
01, --. , %—m are the fundamental roots of A, different from o. Denote the ier 
components of y* and «, respectively by y*, «!”. Then, writing y—1 instead 
of v, from (226) we have 


ee (o euhe A Bie AL Ba i 90), 
v—1 
of!) = of Py TT (0, — 4). 
eK) 
Therefore, putting 
(260) ee = Lf (C1 call 1,) ) 
“x=0 
by (259) we have 
n--m * (LL) (KH) 
(261) ING osu TA 
iil ee 
From (260) 
y—1 
(u) aoe a Pee ssO81% a (1) 
(262) P#\ (6, — 9) =inM! “ ae = 7," 
Since, by (253), the series >) |6,| is obviously convergent, the products 
x=0 
1 0 Ne 
Le tel. o Be 
x=0 Ke 


are also convergent and have values different from zero, since none of the a, — 4, 
vanishes. Further, putting 


(263) I) = ah = I (' se ab 


32” 
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we have (v) raat) ae a) (v+1) 
pe FF =: ie Lae 
u 


co 
Therefore the series }) (Te —- PEED is convergent and O(6,). But this series 


is Tx) 1, and therefore. 
(264) FLO I: 


69. It follows now from (263), (264) and (262) that 


and we have 


(265) N= 3 Pa aatap (1 +00) 
by putting 
d, = sal Caste. 1 ty 
u 
We now put 
JK == Migs fee (hee eo ay 03), 


keep in the sum (265) only the terms for which d,=-0, and combine the terms 
with same value of o,—o. These terms will be now written changing the order 
of the o,, with the index t going from 1 to ¢, and we obtain thus 


NN” t a _ 
(266) Tae aes (AC tO(0, ), 


where now we have added the superscript » to N, obtained from (234) and (261). 
The o, in (266) are now all different one from another. 


We can also assume that 
(267) |o,—o| S|o,—o|S--*S|0,—o]. 
If in (266) we have t=0, then by (254), (256) and (257) we get 
6, O(N EN) = OP OR OOK Pes) 
6, (4 + O(P, rgd?) =e OCR PE a) 
and since P_, K?”-> 0 by (253), 
6, =O (he Pia). 


2.€., (218). Discarding this case, we can now assume / = 1. 


70. Assume now that the condition («) of Theorem 3 is satisfied. Then the 
equality signs in (267) are excluded, and from (266) follows 
{Vo gut daade 
and by (254) and (257) 
(268) Aa Ra 
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and therefore pe 1) 1 o—), 
pil ee ks = oe (x co); 
o—i,, Bp, a-7o 0,—0 


the assertion (218°) follows now with c=. 
Seach 


71, If, finally, the condition («) is not satisfied, but in (267) |o, —o|<|o,—o}, 
then (218°) still holds, and therefore (218°°) holds both for even y and for odd. 
If on the other hand 


|o, —o| =|o,—o 


10) 0 0. = — (0, — 0), 
(266) becomes 


SO ene cna aa, 


ae : (o,—0) (o;—a)") 


But then, since d,;=-0 and d,=-0, either d,+d,=--0 or d,—d,+0. In the first 
case for even vy we have 


Des elasie’s 


Bak (—a)* 


and (218°°) follows immediately. In the second case for odd »v we obtain 


and (218°°) follows again. Theorem 3 is now completely proved. 


72. We apply the method discussed in this part to the matrix (177) treated 
in § 46 by the generalized Rayleigh quotient iteration. Here, using only the 
sequence of vectors €, and introducing again , and o, by (178) and (179), we 
obtain the formulae 

ae Da) Oy ae by Oy 
a) be 26,4 ae, 


Thus, starting with o7=1, we get the following values: 


Oo = 1 Mo = 2 

0, =2 Un AOA? 
Oy = 0.153228 Mz = 405180 
On==,0302723 Ms = .0919831 


04 = .002865057 [a = -007170851 
Os = .0420283226 sts = .0450708. 


For the values of ae y=1,...,5, we have .28, .324, .56, .848, .986. 
y—1 


The convergence to 1 is apparent although not very fast. In order to improve 
the last result, u;, we apply to fs, 4,5 the argument of §19 of Part II, 
replacing the exponents 4 and 3 by 3 and 2 in formula (61). Then we obtain 
the improved value of "5, 45= .0°7998. 
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Erratum 


Berichtigung zu W. Greriscu, Zum Problem der an den Randern fest einge- 
spannten, elastischen Platte, Vol. 2, 227—242 (1958/59). 


Auf S. 240 sind die Fig. 2 und 3 vertauscht. 
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